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The study was concerned with three main areas: 1) analyzing the process of filling a 
fluid management system under orbital conditions, 2 ) determining what portion of the 
filling process requires experimentation, and is subject to modelling, 3 ) developing 
an experimental program plan for evaluating management system concept hardware 
and operational procedures. The fluid management system ^roundruled for this study, 
and identified in Reference 1-1, is a pressure vessel, 1.06m (41.7 inches) in diameter 
with a screen channel device for providing 100 percent liquid flow. Analyses were con- 
ducted or a cryogenic and earth storable management system, where liquid hydrogen 
and N 2 O 4 were selected as representative propellants. 

A propellant fill procedure was designed to meet two objectives. First, system 
pressures must be maintained within acceptable limits without expelling liquid 
propellants, and second, the effective operation of the screen channel device must 
not be impaired by vapor trapped within it during the fill process. 

The case of filling with a cryogenic propellant presents the greatest difficult}- in 
controlling tank pressure without venting liquid. An empty storage tank will reside 
at a substantially higher initial temperature than the cryogen prior to the initial fill, 
and the transfer of energy from the tank to the propellant during the fill process will 
result in venting to avoid exceeding the system structural limits. Because the liquid- 
vapor distribution is not well defined in a low-g condition, venting to relieve over- 
pressure may result in the expulsion of liquid. 
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Vapor will be trapped in the screen channel device in two ways. First, vapor will 
enter the device during the fill process. Second, if helium pressurant is present in 
in the tank, it must be expelled prior to refill, and screen dry out with subsequent 
vapor penetration will occur during the helium vent. 



The filling analysis showed that these problems can be circumvented by introducing 
the processes of tank prechill, fill and tank pressurization to collapse trapped vapor. 

The prechill process is required to reduce tank temperature to an acceptably low level 
prior to initiating the tank fill process. Prechill will consist of a series of charge and 
vent cycles, where either liquid or vapor is introduced during the charge cycle. Vapor 
only will be expelled during the vent cycle because the elevated tank temperature will 
preclude the possibility that liquid will be present at vent initiation. A model, based 
on forced convection heat transfer, was developed to determine the number of charge 
and vent cycles required to achieve prechill. 
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At the completion of prechill, the tank is locked up and liquid introduced through one 
or more spray nozzles to accomplish tank fill. A fill condition of i )0 percent or 
greater will be achieved without the need for venting if near-thermal equilibrium 
conditions are present. It was determined that sufficient bulk fluid agitation will be 
created by the entering liquid to provide near-thermal equilibrium during fill. Together, 
tank prechill and bulk fluid agitation guarantee a no- vent fill. 

Liouid hydrogen system refill is similar to initial fill, once pro-conditioning is com- 
plete. For a system with helium pressurization, pre-conditioning requires that much 
of the helium be expelled before refill can be initiated. The selected procedure 
requires that all usable liquid be transferred to the propellant supply tank, prior to 
initiating vent and purge cycles designed to dilute tank helium mass to an acce blv 
low level. 

Refill will be greatly simplified if helium pressurization is replaced with autogenous 
pressurization. Pre-conditioning is not required in this case. It is necessary only to 
introduce liquid at a sufficiently high flow rate and velocity t^> maintain thermal equili- 
brium during refill. 



At the completion of initial fill or subsequent refills, the tank will be pressurized in 
order to condense the hydrogen vapor trapped within the screen channel device. It is 
estimated that all the trapped vapor will be condensed within approximately five 
minutes, assuring effective operation of the device. 



Fill and refill procedures for an earth storable fluid management system will be less 
complicated than for a liquid hydrogen system. A prechill process will not be required 
because system and propellant temperatures will be approximately the same. Tank 
pressures during fill will remain well below the system vent pressure level. There 
will be no need to employ a purge procedure prior to initiating refill because N0O4 will 
remain within the screen derice during the helium expulsion process. Finally, con- 
densation of N2O4 vapor trapped within the screen channel derice will take an order of 
magnitude less time to accomplish than the condensation of hydrogen vapor. 

The processes selected for further evaluation and, therefore, modelling were prechill, 
fill, and vapor condensation. .All other aspects of low-g fluid management system fill 
and refill were judged to be sufficiently well defined. Preehill and fill are similar in 
one important aspect; it is intended that heat and mass transfer be dominated by forced 
convection in order that these processes remain independent of acceleration environ- 
ment. Consequently, a modelling analysis was performed to identify methods of 
verifying, or modifying, the empirical forced convection equations. This analysis 
showed that forced convection dominance would bo maintained during a one-g preehill 
test and that the results would be directly applicable to a low-g environment. A 
determination was also made that rigorous modelling techniques cannot be employed to 
establish equivalence between normal gravity and low-g fill. It was shown, however, 

xv i 






INTRODUCTION 


The objectives of this study were; to analyze the process of filling a fluid management 
system under orbital conditions; to determine '’•hat portion of the filling process 
requires experimentation, and is subject to modelling; and to develop an experimental 
program plan for evaluating management svstem concept hardware and operational 
procedures. 

1. 1 FLUID MANAGEMENT SYSTEM CONFIGURATION 

Basic to this study is the fluid management system identified in Reference 1-1. This 
system is a liquid hydrogen pressure vessel, 1.06 m (41.7 in) in diameter with a 
screen channel device (Figure 1-1) and helium pressurization system for fluid transfer. 
It contains a high performance thermal control system comprised of a vapor cooled 
shield thermodynamic vent system, multilayer insulation, and vacuum jacket. 

1. 1. 1 SYSTEM THERMAL MASS. The total management system mass of 64.4 kg 
(142.0 lbm) is summarized in Table 1-1. However, an effective thermal mass of 
13.6 kg (30 lbxn) is achieved with the support method of Figure 1-2 which isolates the 
pressure vessel, screen device, and internal hardware from the remainder of the 
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Figure 1-2. Pressure Vessel Girth Joint and Support Method 


system. The thermal mass is defined as that portion of the total which is chilled to 
liquid temperatures during the propellant fill process. For this study, the effective 


mass was assumed to be 13. 6 kg (30 lbm). 

1.2 DESIGN REQUIREMENTS FOR 
CRYOGENIC SYSTEM 

The fluid management system is designed 
for long duration storage of cryogens (7 to 
30 days), during which time fluid is supplied 
to other systems. Depending upon the 
application, flow uemand can be variable, 
as with flow to reaction control system 
thrusters, or well-defined, as might occur 
for continuous flow to a fuel cell. Either 
liquid hydrogen or liquid oxygen will be 
stored in the management system. 

Per agreement with the NASA/LeRC, liquid 
hydrogen was selected as the propellant for 


Table 1-1. System Weight Summary 



k£ 

lb 

— m 

Pressure Vessel 

8.6 

19.0* 

Vapor-Cooled Shield 

4.4 

9.6 

Multilayer Insulation 

4.9 

10.9 

Screen Device 

2.5 

5. 5* 

Collar Supports 

1.0 

2.3 

Internal Hardware./ 

2.5 

5. 5* 

Plumbing 



Vacuum Shell 

19.8 

43.6 

Girth Ring 

20,7 

45.6 

Totals 

64.4 

142.0 


* Mass affected during tank prechill or 
chill. 
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analysis of the cryogenic management system. This selection was made because a 
preliminary evaluation showed that the problems associated with management system 
filling are more severe for liquid hydrogen. Furthermore, the Reference 1-1 design 
was based upon liquid hydrogen as the stored fluid. A system designed for liquid hydrogen 
will also be applicable to other fluids. 

1. 2. 1 LIQUID SUPPLY REQUIREMENTS. Pure liquid must be available at the 

m anagement system outlet at all times. This requirement is translated to mean that 
the screen channel device must be free of vapor once the system is filled and available 
for use. The system configuration (riven in Reference 1-1 was designed without 
consideration bc^g given to orbital propellant fill capability. Consequently, the first 
task undertaken in this study was a selection process to identify those modifications 
needed to fill the channel device with liquid in orbit. 

1.2.2 FLUID STORAGE CONDITIONS. Operating conditions require that liquid 
hydrogen be stored at a vapor pressure of 344. 5 kN/m 2 (50 psta) and that it be 
subcooled by 68.9 kN/m 2 (10 psia) with helium partial pressure. There are applica- 
tions when liquid subcooling is not required. For these cases autogeneous pressurant 
will replace helium. Liquid storage conditions will remain unchanged. 

1.2.3 MAXIMUM PRESSURE ALLOWABLES. Pressure transient will be experienced 
by the pressure vessel propellant tank during the till procedure. Per agreement with 
the NASA/LeRC, maximum allowable pressures cannot exceed 689 kN/m 2 (100 psia) 
at any time during propellant tank fill. 

1.3 DESIGN REQUIREMENTS FOR EARTH STORABLE SYSTEM 

The fluid management system for earth storables will contain cither nitrogen tetroxide 
(NoO^)or monomethyl- hydrazine (MMH). This system will differ from that for 
cryogens in three major aspects; temperature environment, pressurant condition, and 
material compatibility. System temperatures will not vary substantially from propellant 
temperatures at any time prior to or during a fill operation. This is in contrast to (he 
large temperature excursion experienced during cryogen system chilldown. 

System pressures are typically maintained above 300 psia with helium pressurant, which 
will maintain propellants in a subcooled state prior to and following tank fill. There 
are design problems peculiar to the storage of N 0 O 4 which favor titanium rather than 
aluminum as a tank material. But these problems, in general, are related to material 
compatibility and, therefore, are beyond the scope of this contract. Fortunately, the 
thermal mass of the cryogenic system aluminum tank is similar to that of a titanium 
tank that would be designed for earth storables. Because thermal mass is an im- 
portant parameter associated with management system fill, any analysis conducted on 
the aluminum tank will also be applicable to a titanium tank. Consequently, this 
study employed the cryogenic fluid management system for the earth storable fill 
analyses. 
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All analyses are conducted for nitrogen tetroxide (N 2 O 4 ) because its higher vapor 
pressure will make propellant fill more difficult than if monomethylhydrazine (MMH) is 
employed. 

1. 3 . 1 LIQUID SUPPLY REQUIREMENTS. Pure liquid must be available at the 
management outlet at all times. As with the cryogenic system, this means that the 
screen channel device must be free of vapor upon completion of a propellant fill 
process. 

1.3.2 FLUID STORAGE CONDITIONS. The following fluid operating conditions were 
selected for this study phase: 

N 2 O 4 vapor pressure = 130. 9 kN/m 2 (19 psia) 

Helium partial pressure •- 1936 kN/m 2 (281 psia) 

Total tank pressure = 2067 kN/m2 (300 psia) 

1.3.3 MAXIMUM PRESSURE ALLOWABLES. The maximum allowable tank pressure 
for this studj' was selected to be 2412 kN/m“ (350 psiai, which is representative of 
earth storable systems. 

1.4 ANALYSIS APPROACH 

The direction of this study was influenced by attitudes held on manned activity 
in an orbital environment relative to a vehicle system or space station re-supply. 

These attitudes were translated into guidelines to perform filling operations in a 
timely manner and to establish a fill procedure that is relatively independent of 
acceleration environment. 

Regarding the first guideline, that of expediting operations, it is reasonable to assume 
that a fill procedure should be completed within minutes rather than hours because 
system fill is one of many operations to be conducted. However, a reduction in time 
is generally gained at the expense of increased fluid losses, which is undesirable 
due to the h:gh cost of transporting propellants into orbit. Since there was no direct 
method of determining the tradeoff between time and fluid expended, individual 
operations were restricted to approximately 10 minutes, unless the mass expended 
was considered to be excessive. 

Regarding the second guideline, a fill procedure will become more and more independent 
of acceleration if sufficient fluid agitation is provided internal to the tank. That is, the 
heat and mass transfer process to occur during propellant fill will be dominated bv foreei 
convection processes, rather than acceleration, once fluid agitation has increased to the 
level of an inertia-dominated environment. The advantages of an inertia dominated 
environment are clear. First, propellant tank fill will be accomplished in a short time 
period. Second, heat and mass transfer is better defined under orbital conditions if 
internal tank thermodynamics is dominated by the forced convection process. 

Finally, a forced convection environment Increases the likelihood that normal gravity 
tests can be scaled to low gravity. 
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FLUID MANAGEMENT SYSTEM CONCEPT SELECTION 


A comprehensive screening of orbital fluid management system concepts was conducted 
The basic application used was the system described in Reference 2-1, employing 
helium pressurization and screened channels. Systems considered included the L* 

"*** arrangements, pumping, capillary pumping 
using thermodynamm venting, vacuum refilling, inflow baffling, high pressure 
manifolding, shaping of channels and propellant depot refilling stations. These 

*« CT ?r 2 H ,le T or mos ‘ « “or 

in .he following paragraphs are applieable to a 

rangmg from small Spacelab spherical tanks to large vehicle tankage. 


2.1 SELECTION CRITERIA 


f*? ° f tt ! e 00ncepts 3hw ™ * Table 2-1 was compared on the basis of the eleven 
criteria shown in Table 2-2. The comparisons are tabulated in Table 2-3 T h ° 

(Concept 7) to remove vapor from the capillary device. 


2.2 CONCEPT RANKINGS 
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Candidate Concepts for Capillary Device Low-Gravity Refilling 
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Table 2-2. Definition of Comparison Criteria 


Criteria 

Technical 
probability of 
Successful 
Operation 

Adaptability to 

Existing 

Configurations 

Operational 

Simplicity 

Versatility 


Gravity Sensitivity' 
j Ground Testability 


I State of Develop- 
ment/Development 

Cost 

Cost Effectiveness/ 
Recurring 

Weight/Hardware 

Weight/Fluid 


How likely is the concept to accomplish the desired filling of 
capillary devices under all applicable conditions? 

How much modification is required to incorporate the concept 
i„ nr existing capillary device design or configuration? 

operationally, than normal gravity filling o. an empty tank). 

How applicable is the concept to the full range of possible 

eaTllla^ device configurations, including channels, liners, 
start baskets and combinations of these devices? 

How sensitive is the process to low gravity? 
low gravity performance? 

HOW much technology work is required to prove the operation 
of the concept? Will major advances be required? What 
the cost of this development work? 

/ What the cost of building additional units, beyond the 
initial unit? 

What the total weight of all components I device 
modifications required to Implement the system . 

What is the weight penalty due to fluid vented or residual fluid 
that can be attributed to the concept? 


CT J] J j 
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Table 2-3. Rankins: of 


Candidate Concepts 


Key t 2 

Highest - 10 ^ 

Lowest - 0 H 

1 

M 

1 

i 

1 

a 

1 

1 

| 

$ 

!i 

f i 

J a* 

1 

| 

N 

V 

a 
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•» 
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l - 
|S 

flu ■ 

. 91 

n 

* s 

I 3 

111 
3 “ t 
S S i 

t!i 

q je 1! 
J U J= 

0 

3 

3 

N 

1 

* i. 

u 5 
^ 8* 

m > 

1 1 
&J 
1 i 

< (J 

« 

b 

o 

u 

a -5 
£ § 
|1 
i u 

34 
•s ^ 
Ig 

« u 

3 2 

* S* 

a. > 

X3 

o 

l 

& 

< 

t 

! 

O. 

2 

a 

it 
$ i 

1 : 

> s 

i s i 
|i j 
Is* 

C 2 £ 

& t 2 

H U u. 

Open Loop Refrigeration to 
Condense Vapor 

t 

0 

2 

1 
3 
0 

91 

it 

% 

> 

1 

w* 

i 

Shaping Capillary Devices to 
Position and Eject Vapor 

k* 

a 

Ut 

0 

I 

1 

55 

Propellant Control Surfaces/ 
Vent Tube 

3 

Criteria 

2 

3 

4 

9 

6 

9 

7 

12 

10 

11 

1 

* i 

Technical Probability 
of Successful 
Operation 

10 

10 

10 

7 

7 

10 

7 

8 

7 

7 

0 

0 

Adept ability to 
Exiattng Configure- 
tloas 

10 

10 

10 

! 

10 

2 

10 

7 

4 

4 

2 



Operational 

Simplicity 

9 

9 

9 

6 

10 

1 

9 

6 

9 

6 



Versatility 

10 

10 

10 

10 

10 

10 

7 

7 

5 

4 



Oravlty Insensitivity 

9 

7 

7 

6 

10 

7 

3 

6 

4 

1 



Ground Testability 

10 

7 

7 

6 

10 

7 

3 

6 

4 

1 



State of Development 

9 

9 

7 

6 

6 

3 

5 

6 

5 

6 



Cost Effectiveness/ 
Recurring 

10 

9 

9 

S 

5 

9 

8 

6 

7 

6 



Wetght/Hsrdwsre 

9 

3 

7 

9 

2 

l 

7 

4 

5 

5 



Weight/ Vented Fluid 

10 

9 

10 

10 

4 

10 

7 

3 

3 

4 



Reliability /Safety 

10 

10 

9 

9 

10 

7 

9 

8 

10 

10 



Total 

105 

96 

92 

96 

76 

75 

71 

64 

63 

52 

N/A 

N/A j 


demonstrated on the ground. The additional uncertainty with using open loop refrigera- 
tion centers on the need to maintain the liquid in a subcooled condition. Forcing liquid 
into the capillary device through porous tubing (Concept 6) would be straight forward, 
however the difficulties of fabricating u tube of variable porosity caused this concept 
to be downgraded. For the thermodynamic vent/vapor cooled shield (Concept 7) the 
main difficulty would appear to be the possibility of trapping vapor in the channels 
where drag forces would not be able to extract the vapor. The uncertainty with the 
autogenous pressurization concept (Concept 9) would be in controlling mixing, 
pressure reduction and resultant boiling in the capillary device. Use of internal vanes 
to control flow (Concept 10) has not been demonstrated and would require flow analysis 
to determine what tvpes and geometries of barrier would provide uniform flow and 
delay premature wetting. Residuals could also be substantially increased with this 
concept. The use of shaped channels (Concept 11) to eject vapor would impact the 
overall geometry and volumetric efficiency and would, additionally, be sensitive to 
the gravitational acceleration and disturbing accelerations that could adversely 
reposition liquid. Two of the concepts described in Table 2-1 were disqualified 

2-7 ORIGINAL PAGE IS 

OF POOR QUALITY 


(Concept 

ssr '" e (Conoept 31 ■“ - - - “ - Zz?2i77:?°z «• 


2.2.2 ADAPTABILITY OF EXISTING CONFIGTTRATTnwQ ti, j • 

to collapse vapor (Concepts 2 3 4 9 w m * The devices using pressurant 

modifications to "T r6qUire ’ * the mOSt ’ 

added to the eats, tag prcsaurization and fluid acquisition svetems^OnWsT 3 *”!,^ 

= — « 5 L 

revision that will require otabtav .1 T ^f” 10 vent (Concept 7) Is a minor 

liner In probably requiring the addition” *r ° annel devloe but ' vou l d be major for a screen 

entrain trapped vapor Pile tas ld i 1 n , ^ ‘° direCt Uquid flow to °*der «= 

(Concept llTare Tier rewol , ! (C ° nCept 6 > md Sh ^ obannels 

loop refrigeration (Concept 12) would lll^relenslljes' 118 n ° W (Conce P‘ 10 > and °P«" 
configurations. q * extens ive design modifications to existing 


2) and the diffusers (Concept 3) wer^ Hr. 1,11,6 1 pattern - The s Pray nozzle (Concept 

(with the associated plumbine) The th Wng ™ ded because the . v require pressurization 
inlet will require operation otJdto th *™° dynamic vent (Concept 7) vapor cooled shield 
with a dif^^;“““ g C °T ar6d t0 C ° nCept 6 - ThG eonflguratioa 
be switched from the diffuser to the stm S i 00 ^ 1CaCed by the fact that th e inflow must 
to provide a mining £ K 7^^“^ ^ ^ ^ d0Wn; to 0rder 

(Concept 9) is complicated bv ootenM^ * „ Th ? autogenous Pressurization scheme 
mixing after preZrizattan S' IT 7 (°r controlled tank 

control of the disturbances or attitude control” 1 * ° <Condept 111 wtl1 Probably require 
in order to prevent adverse Uqutd/vapor Th^one'T'’ PCrl0dS 

cooling fluid must be used. The proilW dmof l subs ‘antlal quantities of 

more equipment than the other concepts. approac (Concept 8) uses significantly 


vap« 4 tHe C0 " CePl3 p -^^«on to condense 

7 ,he TXT dev,ce USIn * plptas ipside 

vent system flow C ° n ^ rMl0 " S - ™ermody»aniic 

star, baskets. Open loop refr.gera.ton systems ,C«J^^ 1 C 2?Kr i ™ h 38 
devtces without systematical^ reductag iank pressure. 
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concepts (Concepts 10 and 111 have significant influence on the external and Internal 
envelope of the capillary device and therefore may not be directly applicable to some 
liner or start basket configurations* 

2.2.5 GRAVITY INSENSITIVITY. The inflow with spray nobles Concepts 2. 4. ? and 
91 or porous txtblng (Concept di Is Insensitive to^gravity since flow will be inertia 
dominated. The heat transfer between liquid, vapor and tank wall during spraying and 
the condensation of vapor by pressurization will be slightly gravity dependent. The 
flow pattern from a diffuser (Concepts 3, 4 and 91 will also change us a function of 
gravity. Autogenous pressurization (Concept 91 may be additionally sensitive because 
of the Importance of the mixing process. Low gravity heat transfer and tank thermo- 
dynamics are significant factors lit the functioning of the open loop refrigeration system 
(Concept 121. The flow patterns in the internal vane system (Concept 101 used to direct 
flow will lie strongly gravity dependent. Flow patterns, in the concept (concept . i using 
venting to ’’drag" any vapor out of the device will depend upon the balance between drag, 
gravity and surface forces. The shaped channels (Concept 111 depend directly on low 
gravity to operate and will be extremely sensitive to change* in gravity. 

2. 2.(> GROUND TESTABILITY. Ground testability will generally be a direct (Unction 
of gravity sensitivity. Exceptions occur when a process will operate better in low 
gravity than in normal gravity but can still lie demonstrated to operate in normal gra\lt> 
The concept (Concept 2! using spray nozzles Is anticipated to be in this category. 

2.2.7 STATE OF OEYEl.orMEN Y DEVELOPMENT COSTS. The concept using spray 
nozzles (Concept 21 will require a minimum amount of analysis and testing. Diffusers 
and pressurization (Concepts 2. 9. 1. S and 91 add analysis requirement* because of 
the gravity sensitivity of the diffuser flow and the nonequilibrtum pressure rise 
calculations required. Hardware design and fabrication requirements will bo minimal 
for the concepts using helium pressurization (Concepts 2. 1 and 9'. for the other 

concepts, hardware development requirements are more substantial. Development ot 
the porous tube concept (Concept t'-i will tv a difficult fabrication task. No work has 
been done on this concept. Similarly Incorporating shaped channels into a viable design 
(Concept 111 that gives low residuals will be a dttncult design and fabricating task. 

Little work beyond a conceptual nature has been done on this typo of design for other 
than spacecraft applications. Eor the open loop refrigeration system (Concept 121 
development requirement* will include evaluation of heat transfer between the cooling 
fluid, contained liquid and the tank contents and the resulting temperature and pressure 
history. The fabrication requirements of connecting the cc. »irg tubes to the capl\lav> 
device would necessitate a development effort to determine satisfactory attachment 
methods, rhe autogenous pressurization development requirements (Concept 9' are 
greater than that for helium pressurization because of the need for consideration of 
mixing and mixers. Thermodynamic vent system and vapor cooled shield development 
(Concept 7' would require extensive testing to assure that vapor is net trapped. No 
testing of this concept has been done. Internal vanes (Concept 10' would have tv' be 
developed without substantially increasing residuals. Shaped channels (Concept ID of 


1.1 


J JJOAJ.U JJ.'! l.l J l i ITT 


iftAjyiju 



rata 


light enough weight to be competitive would be a serious development difficulty to 
overcome. Developing this concept for typical orbital mission adverse acceleration 
requirements would be a substantial effort. The extensive hardware requirements for 
the propellant depot (Concept 8 ) approach would make this the lowest ranking concept 

for this criteria. 

093 RECURRING COSTS. The recurring costs increase with the amount of hardware 
added to the tank and the complexity of that hardware. The spray nozzle approach 
(Concept 2 and 9) is simplest. The propellant depot approach (Concept 31 is complicated 
by the fittings required to attach the receiver to the depot. Using an inlet diffuser 
(Concepts 3 and 4) is slightly more complicated. Attaching a vent (thermodynamic 
vent to remove fluid. Concept 7) requires additional fittings and valving. Autogenous 
pressurization (Concept 9) may require hardware to control tank mixing after 
pressurization. The hardware requirements are greater for Concept 10) the interna 
vanes inside the channels. Shaped tapered channels (Concept 11) and channels with 
tubing attached (Concept 12) will be more difficult to fabricate. The system with the 
highest recurring cost will be the system using porous tubing (Concept 6 ). The cost 
will increase due to the tubing and the increased structural requirements of beefing 
up the channels to handle the higher pressure and the impingement forces during filling. 

2.2.9 HARDWARE WEIGHT. The spray nozzle concept (Concepts 2 and 9) will have 
the lowest weight, followed by the screened diffuser (Concept 3). The diffuser nozzle 
combination (Concept 4) will be slightly greater in weight. The weight of adding a 
vent to channels to remove fluid ^Concept 7) is next in ranking. Internal vanes (Concept 
10) and shaped tapered channels (Concept 11) will have increased weight because oi 
the added material required to make the vanes and the tapered channel. Tubing weigh: 
(Concept 12) is the main penalty for the open loop refrigeration system. The weight 
penalty for the porous tubing (Concept 0) and channel structural supports penalize the 
system that forces liquid into the channels with variable porosity tubing. Concept 3, 
the propellant depot approach, has significantly more components than the other concepts. 

99 io FLUID WEIGHT PENALTY. In general, minimum fluid weight penalty will be 

due” to the systems using screened diffusers for chilldown and spray nozzles for filling 
(Concept 3) if venting is not reautred during chilldown. If venting will be requires 
during chilldown, the spray nozzle case (Concept 2) gives the lowest vented fluid weight 
penalty. Residuals will not be affected for any of the concepts using pressurization 
('''oncepts 2 3, 4, 3 and 9). For the propellant depot approach (Concept 3) vented 
fluid losses will be a minimum since all vented fluid will be reliquified. The system 
using thermodynamic vent flow (Concept 7) to ’’purge" the channels of vapor will likelv 
have to vent fluid in addition to the normal requirement in order to remove all the 
vapor from the channels. Residuals will not be affected. Using cooling tubes (Concept 
12 ) to condense the vapor in the channels will require substantially more vented fluid 
than the normal boiloff would require. Residuals will not be affected. ResiduaN wi 
be increased for the concepts using porous tubing (Concept 0) shaped channels (Concept 
11) and internal vanes (Concept 10). 
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2. 2. 11 RELIAB ILI TY. Reliability is a direct function of the number of moving parts. 
The systems (Concepts 2, 3, 6, 10 and 11) rated with a score of 10 will have no more 
moving parts than a fluid management system that is not refilled in orbit. They have 
a fill line, pressurization line, outflow line and corresponding valving. The autogenous 
pressurization system (Concept 9) may require systems to control mixing. The concept 
using diffusers and spray nozzles (Concept 4) will have an additional valve (or a three 
way valve) and some controls. The thermodynamic vent system (Concept 7) may have 
an additional valve and controls for operating continuously until filling of the capillary 
device occurs. Open loop refrigeration (Concept 12) will require an additional valve 
and pressure regulator to cool the channels. The orbital propellant depot (Concept 8) 
has considerably more equipment than the other concepts but it is likely to be highly 
reliable and redundant. 

2. £> CONCEPT SELECTION 

Results of the systems comparison indicate that the concepts using pressurization for 
condensing any vapor trapped in the capillary device are the most promising (Concepts 
2, 3, 4and9). However, Concept3will not be given further consideration because 
its primary advantage, that of providing low pressure rise during tank chill, will be 
minimized or nullified due to the introduction of a prechill phase preceding tank chill 
(see Section 3 for further discussion). Concept 4 is also eliminated from further 
consideration because it too includes a diffuser to provide low pressure rise during 
tank chill. Concept 9, which employs autogenous pressurization, is a preferred 
approach because the avoidance of helium can greatly simplify management system 
refill. There are a number of applications, however, which require helium pressuri- 
zation to maintain a continuous supply of subcooled liquid from the fluid management 
system. Consequently, Concept 2 will also be evaluated in greater detail in Section 3. 
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FILLING ANALYSIS FOR SELECTED CONCEPT (CRYOGENS) 


Analyses were conducted for filling the system described in Figure 1-1 with liquid 
hydrogen but modified to allow liquid flow into the tank through one or more spray nozzles, 
and to allow venting through a tapered vent tube. A sketch of the modified fluid manage- 
ment system is given in Figure 3-1. 



Figure 3-1. Modified Fluid Management System With 
Screen Channel Configuration 


The filling of fluid manage- 
ment systems in orbit 
introduces complexities not 
experienced during one-g 
propellant fill. Filling a 
storage tank on the ground 
can be easily accomplished 
because vapor is readily 
expelled as liquid is intro- 
duced. Even cooling a 
storage tank prior to cryogenic 
fill presents no problem. 
Except in the case of complex, 
multiple compartment 
acquisition systems, the 
liquid vapor interface is 
usually well-defined during 
normal gravity filling. Vapor 


will not be trapped within the screen acquisition device if splashing or wicking do not wet 
the screen in advance of the liquid/vapor interface. If vapor is trapped during normal 
gravity fill, its location will be known, and if necessary it can be vented using valving 
and vent lines. 


However, low-g propellant fill cannot be satisfactorily accomplished with standard one-g 
techniques. This is due to the fact that the liquid-vapor distribution is not sufficiently 
well defined to enable pure vapor venting during the filling operation. The influence of 
g-environment upon storage tank fill is illustrated in Figure 3-2. The cost of transporting 
propellants into space for in-orbit propellant fill is sufficiently high that two phase venting 
is undesirable. Vapor trapped within the screen device is also unacceptable because 
pure liquid flow from the fluid management system is required, and cannot be guaranteed 
unless the screen device is free of vapor. ORIGINAL PAGE IS 

OF POOR QUALITY 

Propellant fill is further complicated by two conditions that will exist: an empty storage 
tank will reside at a substantially higher initial temperature than the cryogen prior 
to initial fill, and a partially full tank requires that a major portion of the helium pres- 
surant be vented before refill can be initiated. The question to be resolved for the 
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Fip-re 3-2. Comparison of One-g and Low-g storage 

Tank Pill Employing One-g pm Techniques 
former condition is how tn fin * 

propellants in the process, and without ^xcefd • * ^ Wlth ° Ut ex P endin ? excessive 
pressure. A solution requires that trade* h ^5 structural allowable storage tank 
of tank prechill, chill, and fill. 6 conducted to explore various combinations 

P T o“ storage tank, requires that a ma jo r 

Liquid boiling will occur throughout this helium ^ be taitlated * 

possibility of losing liquid overboard Although tl Whi ° h Wil1 tacrease the 

problem areas, detrirdlnTly^ls^ere longed 16 TV* * the Identified 

system under orbital conditions was devised with / S “ bsequent Procedure to fill the 
variables: liquid inflow rate and velocity i„V T consideration given to the following 
liquid vapor pressure, and liq* d spray cXltfT** temperature > entering 5 
phase of fluid management system fill procedure S CS ‘ DetaIled analyses of each 
sections discussions. are P r °vided in the following sub- 
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3.1 PRECHILL ANALYSIS 

Fluid management prechill is required whenever * 

the stored tank energy will result in * toit al tem Perature is such that 

mod.. PrecbUl [a T"" the •*»— 

taok at velocity that provides good heat^We ° r VaP ° r tot0 the P^P* 11 "* 
walls sod the cooling fluid. This p^ed^Z £ “-Parature 

mass to effeot tank cooling. Liquid hydrogen “ifl be h , reqUirtog U “ le 

than hydrogen vapor in orbit and * ■? T IU ° e readll Y available 

effecting prech,/ ,t is ,^sWe how,ver t r. h ‘° ^ “* fla “ CandldMe »* 
the transfer lines to the pressure vessel cL l VaP ° r nonnall y generated in chilling 
is given in this section on prechiU "“p^ ** PreChiU - A ^—sion 

3.1.1 PRECHILL WITH HYDROGEN VAPOR Thp 

prechill is to reduce tank temperatures' s^fficien^v th requirement fo ** system 

accomplished without venting. K is illicit chiH “d will be 

unacceptable during the chill and fill ? u ^ requirement that venting is 
unknown quantity of liquid will be lost o^rto^stace f the n P ° Ssibilit 5 r *** an 
maintained during this process. Venting ST COntro1 cannot be 

because ihe elevaled lauk ,mpera.u re ^ ^SS^Xphase. 

ments is th^orS m^um pressurTSaf to determinto e PrecWll require- 

functiou o f initial tank temperature Maximu^ T °° 0Ur durt " ? ta ” k chm as a 

addition, occurs when the tank vapor anTtTJT 1? pressure > for a given mass 

The analytical development for theoretical maurim reSide “ tke same temperature. 

tneoretical maximum pressure is given below. 

The Firs. Law expression for introducing l,,„, d , Dt0 . contalner y 
dEg + dE w = h^drtiL 


dEg = (U g m g , 2 - ( Wl , change id ullage energy 
If one assumes that the tank is initially evacuated 


(3-1) 


m g 1 = 0 and dEg « u. 


«2 m 82 


Also for an initially evacuated container, m g2 . dm L 

dEw (“ivmw>2 - (UwUMj = change in tank wall energy 
Since tank mass is constant, dE w . <Uw2 - u^, n^ 

Combining Equations 3-1 through 3-4 
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Finally 


U ?2 m g2 + (Uw 2 ~ “wi> “w = h L ^g 2 
(U S2 " h L> m g 2 = - «w 2 ) “w 

m S2 = (U w x - V 2 > m ^v/(Ug 2 - h L ) 

where 

“w 2 and u g 2 are evaluated at T 2 

dE g = change in tank wall internal energy 

h L = enthalpy of liquid entering tank 

dm L = differential liquid mass addition to tank 

u g = internal energy of vapor in tank 

mg = mass of vapor in tank 

°w = Internal energy of tank wall 

“w = tank wall mass 

T = temperature 

subscript 

1 « conditions at beginning of interval 

2 = conditions at end of interval 
From the equation of state, gas pressure is 



/here 

2 = compressibility factor 

R = gas constant 
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s tank volume 


- gas pressure 


„ , j flformine d as a function of mass addition by a 

Peak pressure during task cWU«.l» _ Results are * Flgu re 3-3 for a 

simultaneous solution of Equ initial tank temperature conditions. Note 

range of incoming liquid vapor pressure and initial tan ^ gures of this figure 

.. . , Hwt na 


(120) 1 
800l 


INITIAL tank TEMPERATURE - 
. 194. 4K (350R) 


138. 9K (250R> 


mat tuc r- 

are the theoretical maxim urns that can 
occur for the given initial tank 
temperature. A lesser or greater 
mass addition will result in a lower 
tank pressure than the theoretical 

maximum level. Note also that 

entering liquid vapor pressure will 
have a minor influence upon tank 
pressure. 


y::tering liquid 

VAPOR PRESSURE* 
kN/m- (P«U1 — 


■ . 275 . 6 ( 40 ) 

137.8 ( 20 ) 


( 2 ) ' ’ 
HYDROGEN MASS ADDITION. Wg i^W 


Figure 3-3. 


Peak Pressure During Manage- 
ment System Chllldown Versus 
Mass Addition and Entering 
Liquid Vapor Pressure 


The theoretical maximum pressures 
described by Figure 3-3 are plotted 
in Figure 3-4 as a function of initial 
tank temperature. Initial temperatures 
in excess of 194. 4K (350R) were not 
considered because of the requiremen 
t to maintain system pressure below tl* 
maximum allowable level of 689 kN/m- 
(100 psia). This figure shows that 
maximum pressures will be 50 percent 
or less of the maximum allowable for 
initial tank temperatures lower than 
138. 9K (250R). 


. .wiling the oressure vessel to the lowest temperature 
There le an obvious advantage to chill over-pressure can be eliminated, 
possible because any concern rel ^ complication Incurred as more hydrogen 

This advantage must be weighed g sy3Um t0 lo wer and lower temperatures. An 

vapor and time Is required to pre ' , c hilllng to liquid hydrogen tempera- 

extreme condition to consider Is requirements will be quite 

lures. Without analysis one can , tempera ,ures and, therefore, vent 

higi*. at low tank temperatures beca ^*v-. rmore the likelihood of liquid residing in 

enthalpy will be correspondingly ow- U ’ which increases the probability 


u«n ivrT rioi A Drechill subroutine was Introduced into the HltPRE 
3. 1. 1. 2 Tank Prechill Model the varlab ies of a prechill procedure, 

computer program in order to adequately as 
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me subroutine is described 
in greater detail in Appendix 
B, whereas the model 
assumptions, results and 
conclusions are given in 
the following paragraphs. 


(* 00 ) 


(MO) 


(300) 


INITIAL TANK TEMPCRATURE, K (R) 


Figure 3-4. Maximum Pressure During Management 
System Chilldown Versus Initial Tank Temperature 
and Entering Liquid Hydrogen Vapor Pressure 


The key to this analysis is 
in identifying the appropriate 
heat transfer coefficient 
between incoming vapor 
and the tank walls. 
Conventional forced 
convection expressions 
for flow over a flat plate 
or for jet impingement 
(MO) u P° n a surface, represent 
configurations that are too 
dissimilar to be applicable. 
Instead, a correlation 
developed for industrial 
mixing processes was 
selected as being 


required modifications are given as ' ' ™ h * tnuisfer correlation and the 


JL L V 3 -n f ( p / V)M V /4 

PC P V PR / * 13 1 p 2 I (Reference 3-1) 

vhere 

h = heat transfer coefficient 

o " fluid density 

c p = constant pressure heat capacity 

n PR - Prandtl number 

P[ = mixer input power 

V 


M 


= tank volume 
= fluid viscositv 



Equation 3-9 was developed for liquids contained In cylinders. These liquids were 
continuously agitated with a mixing unit. Mixer input power was responsible for fluid 
agitation and is one of the variables of Equation 3-9. For the prechill application , it 
is believed that fluid agitation will be the same whether a mixer or fluid inflow is 
responsible. It is important only to have equivalent power conditions. Power output 
rather than power input will influence fluid agitation. Consequently equivalence will 
be between fluid power input and mixer power output. This results in 


P t (EFF) = P Q = mr 


(3-10) 


where 


P Q = mixer power output 

EFF = mixer efficiency (conservatively assumed as 40 percent for this study) 
m = entering mass flow rate 

v = entering fluid velocity 

mv 2 = fluid power input 

Substituting Equation 3-10 into 3-9 results in 


h 

pC P 


(%r) 


2/3 


0. 163 


(m v 2 /V) 



(3-11) 


Equation 3-11 indicates that heat transfer to the tank walls can be controlled by varying 
entering flowrate and velocity. 


The following charge and vent procedure was selected for fluid management system 

prechill 

1. Charge the tank at a known vapor flowrate until the difference between wall and gas 
temperature has reached a specified value. 

2. Vent the tank to a pre-detennlned low pressure. A reasonable level was selected 
as 6. 89 kN/m 2 (1. 0 psia). 

3. Charge and vent the tank as required to reduce tank temperature to the pre-determined 
level. 


Figure 3-5 gives prechlll mass as a function of final tank temperature. Note that prechill 
mass requirements Increase as final tank temperature decreases. Also, prechlll mas 
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Figure 3-5. Hydrogen Vapor Mass Required to Prechill 
the Fluid Management System Tank 


can be reduced by decreas- 
ing charge flow rate or 
increasing vapor velocity. 

As expected increased 
prechill times are required 
to provide lower final 
temperatures, Figure 3-6. 
These times can be reduced 
by increasing charge flow 
rate or by decreasing vapor 
velocity. Finally, Figure 
3-7 shows that the frequency 
of charge and vent cycles 
increases as final tank 
temperature decreases. 

This figure also shows that 
fewer cycles are needed 
for tank prechill as charge 
flow rate is increased. 


a 

| 350 

y 

ui 

Vi 


u 

u 


250 


L50 


An evaluation of Figures 3-4 through 3-7 indicates that a final tank temperature of 
about 138. 9K (250R) represents a reasonable trade between the advantage of low peak 

pressures during tank 
chill and the disadvantages 

450 L*"' — tSr * j c-iTic.iui’iu on;- 1 c.mr. - oo.or*. jovn». | of increased hydrogen 

mass, time, and charge 
and vent cycles. These 
figures also indicate 
that (1) the influence of 
entering velocity is 
minimal, (2) prechill 
mass requirements are 
directly proportional to 
entering flowrate, and (3) 
prechill time and charge 
and vent cycles are 
inversely proportional to 
entering flowrate. 
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System With Hydrogen Vapor 

3-8 


3. 1.1. 3 Summary. The 
above parametric data has 
resulted in the following 
recommendation for system 
prechill: 
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TIMli FROM START OF PRECHILL (SECONDS) 


Figure 3-7. Fluid Management System Tank Wall Temperature 
History During a Gaseous Hydrogen Prechill 


prechill flowrate = 0. 0045 kg/sec (0. 01 lb m /sec) 
incoming velocity = 61 m/sec (200 ft/sec) 

The tank will be prechilled to 138. 9K (250R) and 0. 66 kg (1. 45 U^) will be expended. 

GHo will be expended during a 295 second period. The higher flowrate was selected 
because it will reduce the required number of charge and vent cycles from nine, at. a 
flowrate of 0.0023 kg/sec (0. 005 ll^/sec), to five. The increased prechill mass resulting 
from the higher flowrate is acceptable, considering that the smaller number of charge 
and vent cycles will simplify the procedure. 

3. 1. 2 PRECHILL WITH LIQUID HYDROGEN. It is expected that liquid hydrogen will 
be more readily available in orbit than hydrogen vapor. Thus it will be more convenient 
to initiate prechill with liquid. Thermodynamically and fluid dynamically there is little 
difference between a liquid and vapor prechill. Approximately 15 percent less hydrogen 
is needed for prechill when it enters as a liquid, due to the heat of evaporation which 
increases its energy absorbing capability. Heat transfer coefficients at a given fluid 
power input, will be the same as for gas phase flow, once the liquid is evaporated. The 
primary difference between liquid Inflow and vapor inflow results from film boiling 
which occurs when liquid strikes the hot tank walls. But, since film boiling heat 
transfer rates will represent only approximately 15 percent of the fluid energy absorbing 
capability, total prechill time and the number of charge and vent cycles will not differ 
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significantly between liquid and vapor inflow. It is concluded, therefore, that the 
results for prechill with vapor inflow apply equally to a liquid inflow. 

3.2 TANK CHILL ANALYSIS 

The fluid management system tank temperature will not exceed 138. 9K (250R) at the 
start of the chill process. Liquid hydrogen flow to the tank will be initiated once the 
system vent valve is closed. Film boiling will occur as liquid strikes the warm tank 
walls, and the resulting evaporation will increase absolute maximum pressure to 339 
kN/m 2 (49. 2 psia), as indicated by Figure 3-4. This pressure will be about 30 percent 
of the maximum allowable level. Neither tank chill nor the subsequent tank fill process 
will be influenced whether pressure peaks at 339 kN/m 2 (49. 2 psia), or at a lower level. 
Asa result, it will not be necessary to impose constraints on the flow process for 
purposes of controlling tank pressure and, it will be possible to perform tank chill 
with any constraints imposed by either prechill or tank fill.. Thus it is seen that the 
selected piechill process has virtually eliminated all concen s relative to an in-orbit 
tank chill process. 

3.2.1 TANK CHILL THERMODYNAMICS. Despite the knowledge that a no-vent tank 
chill can readily be accomplished, it is still necessary to evaluate the influence of 
inflow parameters upon tank pressure during this process. Tank pressure determination 
Involves a complex interrelationship between the incoming liquid, ullage mass and tank 
mass. A computer program (HYPRES), which is described in Appendix A, has been 
written to analyze the tank chill process. Briefly, it is assumed that liquid enters the 
tank as a spray, and exchanges energy with vapor as it flows through the ullage. Heat 
exchange between the wall and ullage vapor is also permitted. Liquid droplet evapora- 
tion or condensation will occur, depending upon liquid spray characteristics and 
ullage-to-liquid temperature differences. Liquid-to-wall heat transfer is limited to 
film boiling. The nucleate boiling process is not included because the stored wall 
energy is insignificant in this temperature range. 

Two variables requiring further analysis are the drop size determination and wall 
boiling during the chill process. Droplet size determination is required to verify that 
spray conditions reside within the acceptable analysis range. Wall boiling rates must 
be Identified as a function of such variables as fluid properties, velocity, drop size and 
gravity* 

3.2. 1. 1 Liquid Droplet Size Determination. Spray-nozzle drop size information is 
available from manufacturer's data for air-water combination, but rarely for other 
fluids. No all-purpose defining equation for the effect of physical properties exists. 

The dependence of drop size on geometry, velocity, and properties is complicated by 
changes in the nature of the break-up process as these variables change. The 
suggested predictive route starts with manufacturer's air-water data at the desired 
flow rate, and then drop size corrections are made for different physical properties. 

The recommended equation for drop correction is given in Reference 3-2 as 




zmn :~i ~i ~~j 


.1L-L 
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dvm 
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dvm r c 

n 0.5 r M 
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(3-12) 


where 

dvm volume median droplet diameter 
a => surface *ension, dynes/cm 

Mi = liquid viscosity, centipoise 

= liquid density, lb^/'ft^ 

The exponential dependences in Equation 3-12 represent averages of values reported 
by a cumber of studies. 

A more significant drop size than dvm is the Sauter mean diameter, d, which has the 
same ratio of surface area to mass as the total drop population. This diameter was 
used in the study because a representative surface area is critical to the chill process. 
It is typically 70 to 90 percent of dvm; a value of 80% was selected for this study. 

3.2. 1.2 Spray Nozzle Selection. A full cone spray nozzle was selected for 
propellant tank chill and fill. Typical nozzle performance characteristics were 
obtained from Spraying Systems Co. , Including particle size as a function of nozzle 
inlet pressure for a variety of nozzles, Figure 3-8. Equation 3-12 was employed to 
determine drop size for the combination of liquid hydrogen flowing through gaseous 
hydrogen. Droplet diameter is plotted versus liquid flowrate for two spray nozzle 
configurations, Figure 3-9. These configurations were selected because of the large 
drop size c- ■sated during the spray process. Note that the maximum hydrogen drop 
size is expected to be less than 600 ju (0. 024 inches). 


3.2. 1. 3 Liquid Splattering on a Hot Surface. The actual heat exchange process to 
occur during tank chill will be that of heat transfer between an individual droplet and the 
hot tank walls as contact occurs. This method of heat exchange is different from that 
of pool film boiling for one obvious physical reason; film boiling is based upon a hot 
surface Immersed in a pool of liquid where a stable vapor film has developed, whereas 
heat transfer due to droplet splattering is of a more transient nature where a vapor 
film is created and destroyed with each impact. The following expression was 
developed in Reference 3-3 for the maximum heat transfer experienced by the Impinge- 
ment of a single droplet upon a hot surface 


^max 

P L d 3 X 


8. 44 x 10“ 4 



0.341 


(3-13) 


and 
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Figure 3-8. Water Spray Particle 
Size Versus Pressure For Fulljet 
Nozzles 


’I’f = (T w +T s )/ 2, film temperature 

o' = surface tension 


Qmax = maximum heat transfer per drop 

P L = liquid density 

d = drop diameter 

X = modified heat of evaporation by 

. Equation 3-14 

hfg = enthalpy of evaporation 

Cp y ~ vapor specific heat at constant 
pressure 

T w = wall temperature 

T s = liquid saturation temperature 

V = drop velocity 

PV F = vapor density evaluated at T p 


Equation 3-13 was based upon results of tests with water, Acetone, and Ethanol. It was 
observed that the maximum heat transfer condition corresponded to a droplet impact 
angle of 27 degrees as measured from the normal to the hot wall. The equation is valid 
for values of the independent variable, PL 2 V 2 d/p VF J c , less than loL Furthermore, 
it was determined that the efficiency of the process, defined as the ratio of actual heat 
transfer per drop to total energy required to evaporate the drop, is of the order of 50 
percent for a single collision. Subsequent secondary collisions of splattered particles 
with the hot surface tend to increase the overall efficiency. 

Observations made during the Reference 3-3 experimental effort indicated that the heat 
transfer mechanism was that of conduction across a vapor film. Qualitatively, a droplet 
will begin spreading on impact with the hot surface and the combination of contact 
surface area. A, and film thickness, 5 , with time influences heat transfer rate. Both 
A and 5 are influenced by drop size, impact velocity and impact angle, although the 
exact relation is not known. 
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UOUID SPRAY FLOWRATE, lcg/»«c <W»ec) 

Figure 3-9. Liquid Hydrogen Droplet Diameter Versus Flowrate 
for Two Spray Nozzle Configurations 

Maximum heating rate for a spray condition can be determined by combining Equation 
3-13 with the following expressions 


Q_ n „ = n x Q 
max ^ 


max 


(3-15) 


m = n x (mass/drop) - n 
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P L d3, 
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m x Q 
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^max ~ n 3 - im\) x 8.44 x 10 (« 2 V d/oy z) " (3-17) 

P d -Li F 

6 m L 


0.341 


= 4.42 x 10 " 4 X (p L 2 V 2 d/Py r ) 


0.341 


(3— IS) 


* 

^max = maximum heat transfer rate under spray conditions 
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n * number of droplets per second striking a hot surface 
m - mass flowrate of spray condition 

Fquntton 3-18 luis been Introduced Into the HYPWKS computer program . The significance 
of this expression Is that the splattering phenomenon appears to be Independent of 
acceleration environment. Thus It is reasonable to expect that chill tests, conducted 
In a normal gravity environment, will bo applicable to a near 7 , 0 m gravity environment. 
This is not to suggest that scaling or extrapolation problems mav not exist. However, 
utv ertnlntles In wall boiling heat transfer rate should not invalidate the liquid spray 
approach to tank chill. Furthermore, regardless of the uncertainties of Kquatlon 3-18 , 
peak pressure will not exceed 335) kN/m- (49. '2 psla). 


3.2. 1.4 Peak Tank Chill Pressures. Computer runs were conducted to obtain para- 
metric data on the Influence of liquid flow variables upon maximum pressure during 
tank chill. Figure 3-10 provides a summary of these computer output results. The 
significant results of this Figure are (U peak pressures Increase in magnitude at 
’ 'dueod Inflow rates, and (21 peak pressures are Insensitive to flowrate, droplet 
velocity and droplet diameter at flowrates exceeding 0. 18 kg/ see (0.40 ll>j n sec). These 
observations can bo explained as follows: maximum pressure during tank chill occurs 
when all the tank wall energy Is eonveeted to the ullage vapor. The other pressure 
extreme occurs when tank wall energy Is absorbed by liquid. Thus, peak pressures 
occur at low flowrates because much of the liquid Is evaporated before striking the 
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tank wall with the result that considerable wall-to-ullage energy transfer ooours. 
Conversely, the high flowrate oondltlon will mulntaln a liquid film on the tank walls 
and absorb all the wall energy. 

3.2. 1,5 Summary. Analysis shows that the selected prechlll prooedure has virtually 
eliminated all conoem relative to a tank ohlll procedure. This Is because the theoretical 
maximum pressure to oocur during a no-vent tank ohlll will be about 3,19 kN/m - (49. - 
psla). Figure 3-10 Indicates that the expected maximum pressure will be less than 27« 
kN/m 2 (40 psla) If liquid Is sprayed Into the tank. There Is no possibility that venting 
will be required since the maximum allowable pressure Is (589 kN m- (100 petal. 
Furthermore, any uncertainty In film boiling heat transfer rate that may exist due to 
liquid Impingement on a hot wall, cannot possibly result lit a pressure greater than the 
theoretical maximum quoted above. 


3.3 TANK FILL ANALYSIS 


Propellant tank fill represents a continuation of liquid flow Initiated for the tank ohlll 
prooess. Tank fill Is defined as commencing at the Instant that tank temperature drops 
to liquid temperature. This definition has been selected for the convenience of analysis 
since from this time on the fill process Is not Influenced by fluid managenwnt system 
configuration, other than volume. That Is, propellant state conditions during fill are 
only a function of the thermodynamic Interaction between liquid and vapor. 

3.3. 1 THERMAL EQUILIBRIUM TANK FILL. The single requirement for tank fill Is 
to maintain acceptably low pressures during the process. Tank pressures will be at a 
minimum If thermal equilibrium conditions are maintained during fill. Thermal 
equilibrium will be approached as heat and mass exchange between the phases Is 
increased, which can be achieved bv creating a highly agitated fluid condition. 

Given the assumption of a thermal equilibrium tank fill one obtains a simple relation- 
ship between Initial tank temperature at the start of no-vent chill and fill. Incoming 
liquid vapor pressure, and final liquid vapor pressure. This relationship Is derived 
below from the First Law of Thermodynamics for liquid flow Into a -closed container 

dKg ♦ dE l, • dK w • dmj (3-191 

dF.r a tujmA* - pumA, change In ullage energy (3-201 

d K i (Upm^)-j - (Uj raj -=> change In liquid energy (3-211 



hj dmj •> hu m L» ‘ m l enthalpy change due to entering liquid (3-221 
dh\ v - (W 2 ~ («\v m \vh n 'w i u 2 -u l\v 


» m w '.u^y change In tank wall energy 
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If one assumes that the tank is initially evacuated 


m gl - “L, * 0 


and 


dE g ' “*2 ”g 2 


dE L = u L2 m L2 
h L dm L = h L m L2 

Combining Equations 3-19 through 3-26 we have 


a g 2 m g 2 + u L 2 “L 2 * "V ‘ “w - h L “L 2 
Solving Equation 3-27 for hL, results in 


h L ’ °ft m g 2 /m ^2 + “L 2 * lu w "V'“L 2 





At thermal equilibrium one can relate hL» uj^ and Ug^ to liquid vapor pressure and 
temperature. Consequently, UL 2 and Ugg are known once final vapor pressure is 
specified. can also be determined since initial and final temperatures are given. 
Finally, hL (which is a function of entering liquid vapor pressure) can be determined 
for a desired liquid fill condition. 


Equation 3-28 is summarized in Figure 3-11 which gives entering liquid vapor 
pressure as a function of initial tank temperature and final tanked liquid vapor pressure 
for a 90 percent and 100 percent liquid fill condition. Note that final liquid fill 
condition has a minimal influence on the other variables. Note also that final vapor 
pressure is less than entering liquid vapor pressure. This difference is due to the 
combination of initial tank wall energy and the heat of compression which are released 
to the tank fluid during chill and fill. 


3. 3. 2 FILL MODEL FOR LIQUID SPRAY DOMINANCE. The intent of the tank fill 
process will be to create conditions conducive to attaining near-thermal equilibrium. 
These conditions will be achieved by introducing liquid into the tank through spray 
nozzles, which resulting spray will create a large liquid-vapor surface area. The 
combination of large surface area and high droplet velocity will provide the high heat 
transfer rates needed to attain near-thermal equilibrium conditions. As tank fill 
continues, the internal tank fluid environment will change from one where liquid droplets 
reside within the ullage volume to one where vapor bubbles are entrained within a 
liquid bulk. The transition from heat transfer dominated by liquid droplets to heat 
transfer dominated by vapor bubbles cannot be determined precisely, but it is 
expected to occur in the range of 10 percent to 20 percent liquid fill. 
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TEMPERATURE AT START OF TANK CHILL, K ;R) 

Figure 3-11. Fluid Management System Final Tank Pressures 
for a Thermodynamic Equilibrium Chill and Fill 
Process 

3.3. 2.1 Initial Fill Pressures. Equations describing heat and mass exchange 
with the ullage are those given in Section 3.2, and contained in HYPRES.for the tank chill 
process. 

Propellant tank pressure and hydrogen temperature histories have been determined using 
the HYPRES computer program and results are given in Figure 3-12 and 3-13 for inflow 
rates of 0.091 kg/sec (0.2 lbm/sec) and 0.227 kg/sec (0.5 Ibm/sec), respectively. Note 
that there is a gradual pressure decay rate once peak pressure has occurred. This Is 
due to the influence of incoming liquid as it continues to chill the ullage. Also, it is 
seen that vapor temperature approaches liquid temperature within seconds from the 
start of liquid Inflow. The temperature histories show that near-thermal equilibrium 
conditions are achieved early in the tank fill process. 

3.3.3 FILI. MODEL FOR VAPOR BUBBLE DOMINANCE. As the propellant tank begins 
to fill with liquid, the mechanism of liquid spray in a vapor environment will change 
to that of vapor bubble entrainment and dispersal within the liquid volume. Entrainment 
will occur as a result of liquid impingement with the liquid bulk. Dispersal will be 
caused by agitation created by liquid sprayed into the tank. It is estimated that this 
mechanism will be the dominant mode of heat and mass exchange throughout most of 
the tank fill process . 
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SYSTEM PARAMETERS AND FLOW CONDITIONS 

1. TANK VOLUME * 0. 62 M 3 (22 FT 3 ) 

2. TANK MASS » 13.6 KG (30 LDm) ALUMINUM 

3. INITIAL TANK TEMPERATURE * 138. 9 K (2S0R) 

4. ENTERING LIQUID VAPOR PRESSURE » 137.8 kN/M 2 (20 PSIA) 

5. LIQUID SPRAY FLOWRATE - 0. 091 KG/SEC (0. 2 LBm/SEC) 

6. DROPLET DIAMETER » 762 u (°- 06 DO 

7. - DROPLET VELOCITY - 24.38 M/SEC (80 FPS, 
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Figure 3-12. Fluid Management System Pressure and Temperature Histories During 
Initial Ten Percent Fill Period (Low Flowrate) 


SYSTEM PARAMETERS AND FLOW CONDITIONS 

1. TANK VOLUME ’ 0.62 M 3 (22 FT 3 ) 

2. TANK MASS * 13.6 KG (30 LB M ) ALUMINUM 

3. INITIAL TANK TEMPERATURE = 138.9K (2S0R) 

4. ENTERING LIQUID VAPon PRESSURE - 137.8 kN/M 2 (20 PSIA) 

5. LIQUID SPRAY FLOWRATE » 0.227 KG/SEC (0.5 LBm/SEC) 

6. DROPLET DIAMETER* 762 (i (0.06 IN) 

7. DROPLET VELOCITY - 24.38 M/SEC (80 FPS) 
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3. 3.3. 1 Vapor Dispersal Mechanism. There are many Industrial applications for 
which vapor dispersal in liquid Is essential to achieving a high rate of heat and mass 
transfer. This subject has been extensively studied and much of the work has been 
collected and summarised in Reference 3-1. The dispersal of vapor in liquid is 
brought about by fluid dynamic forces which have to overcome the static forces of 
surface tension. Such surface forces resist dispersion by attempting to retain bubble 
sphericity and prevent gross distortion leading to break-up. The dynamic forces 
which bring about dispersion may be due to bouyancy or induced fluid flow creating 
viscous or inertial force which, if they do not act equally over the surface of a bubble, 
may cause it to deform and eventually break-up. Mechanical agitation is employed 
in mixing vessels to create shear stress by means of turbulence, which in turn causes 
bubble break-up. The empirical equations describing bubble diameter, percent vapor 
volume entrained beneath the liquid surface (vapor hold-up) and heat and mass transfer 
coefficients are applicable to a configuration similar to that shown in Figure 3-14. The 
figure describes a cylindrical tank filled with liquid, within which is immersed a 
mechanical device for mixing the liquid. Vapor is introduced into the tank through a 
manifold located near the bottom of the tank. The large vapor bubbles that form near 
the manifold are broken up into numerous small bubbles as a result of liquid shearing 
forces counteracting and overcoming surface forces of the individual bubbles. This 
bubble breakup mechanism is a Weber number phenomenon. 


The equation for determining bubble diameter is given in Reference 3-2, Equation 
13-138 as 


a no™ °' 6 °* 5 

0. 0279 o j 


0.2 , /,, x 0 - 4 

V (Pi/Vr) 


+ 0.09 


(3-29) 
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Figure 3-14. Apparatus for Study of Vopor 
Dispersal in Liquid 


3-19 






vi .i- 1 ! ) j. j j j j : jj j j j ~] ~ttt 




r* 



P L = liquid density 
Pjj20 = densit y °f water 

Pi = input horsepower to the mixing device, hp 
By expressing p in watts and V L in »*. bubble diameter can be expressed as 


d » 


1.636 a 0 * 6 c °- 5 


0. 2 . 0. 4 

y (Pi/v L) 


+ .09 


(3-30 


thosf due^to^gi^ equated "* 

described below, must be aonliprt tm quated Wlth the m so a separate treatment, 

those which determine heat L ma Jf T ’ which treats turbulent forces as 

bodies submer^r ta mLta vessefo 1Z ^ aISO »*“ Wiled to fixed 


to a d,spersiOT —« * 

may be used to give Mo~ oolhe t b . “ “* flU ‘ d ' The the ^ •» local isot, 
around the partfole wtcT?esuU s ,l th “, ,f ta the sma11 fluld volume 

number to the foUowta S o^ession for turbulence Reynold 


N Re * » L 1/3 (P|/V L ) 1/6 d 2/3 /p 


1/2 


where 


(3-31) 
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where Npr a Prandtl number, and units for the variables h/(CppL) anc * 
rr/n./VT )Ul3/pt 2 1 are g* ven to m / sec< Equation 3-32 Is applicable over the 

range 10-4 to 10 2 for the parameter {[(pl/VL) Hl3/°L 2 3 1/4 1 N p r 2/ * Thls 
equation is Identical to Equation 3-9 except that liquid Is used. 

3. 3. 3. 3 Bubble Diameter During Tank Fill. Equation 3-29 cannot be directly employed 
to determine bubble diameter during the tank fill process because the expression 
includes power Input to a mixing device which agitates the fluid. However, substitution 
of Equation 3-10 Into 3-29. and once again assuming that mixer efficiency is 40 percent, 

results In 


, °- 6 °- 5 
1. 134 a ( 

0.2 . 2 , . 0. 4 

v (mv /VjJ 


+ 0. 09 


(3-33) 


Bubble diameter is plotted versus € and mv 2 in Figure 3-15 for saturated liquid 
hydrogen properties evaluated at 276 kN/m2 (40 psla). The decrease in bubble 
diameter as mv 2 Increases Is expected because the Increased turbulent forces will be 
balanced by the higher surface tension forces of smaller bubbles. Bubble diameter 
will Increase as € Increases due to the Increased rate of bubble coalescence as the 
volume percent of entrained vapor is increased. 

3. 3. 3. 4 Vapor Entrainment During Tank Fill. Vapor holdup, f must be known before 
Equation 3-33 can be employed to determine bubble diameter. For a typical industrial 
application, the vapor inflow rate Is an input variable. During tank fill, however, vapor 
entrainment will be caused by the interaction of liquid jet with liquid bulk as one 
penetrates the other. The following simplifying assumptions have been made in order 
to determine the quantity of entrained vapor In liquid: 

1. Acceleration environment results in a nearly flat liquid-vapor interface. 

2. Mixing between incoming fluid and liquid bulk is for a circular jet. 

3. Mixing between re-entering jet and liquid bulk is for a circular jet. This allows 
for a rapid jet velocity decay and, therefore, decreased vapor penetration depth. 

4. Entrained vapor flows at the average jet velocity, v b x* 

5. Vapor penetrates liquid to a depth, X, at which time v LX = 10 v b , where v b - bubble 
rise velocity. 

6. Bubble stay-time, t = x/vt>. 

Figure 3-16 illustrates a mechanism for entraining vapor during the fill process. A 
liquid jet is Introduced into the tank and entrains additional liquid as it flows through 
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Figure 3-15. Sauter Mean Diameter of Hydrogen Vapor Dispersed 
In I iquld Hydrogen Versus Input Fluid Power 
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Figure 3-1 fi. Meclianlsm for Entraining Vapor In Liquid 
During Tank Fill 
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the bulk toward the ullage. The jet, which has Increased in size, rises about the liquid 
surface, traverses the ullage space, strikes the opposite end of the tank and is deflected 
back to the liquid. Vapor is entrained into the liquid bulk as a result of jet penetration. 


Vapor entrainment in liquid can be estimated as 


w 


(3-34) 


The volume rate of vapor entrainment can be approximated from the following empirical 
expression (Equation 18-133, Reference 3-2) 


fp L\0.5 f 

Vg= 0-0316 [jT-J ^ L *v l 


i, * V 


= 7.1VI 

(for LH 2 properties evaluated at 276 kN/m 2 (40 psia) 
where 

Vg -- volume rate of vapor entrainment into liquid pool 
L = length of liquid jet (defined by Figure 3-16) 

Vl = volume flowrate of jet re-entering bulk 
V L = velocity of jet re-entering bulk 

a = surface tension 


(3-35) 

(3-36) 


It is reasonable to assume that Vl and v L remain constant as the jet travels through 
the ullage before returning. It Is now necessary to define V L and v L as a function of 
the entering jet properties and tank geometry. 


Jet entrainment can be described as a constant momentum process. Thus mv (entering 

jet momentum) = m L vl (exiting jet momentum). Furthermore, mv = p V *v = a v, *v T 
Or L L - 


V L x v l = mv/P 


Substituting Equr.tion 3-37 into Equation 3-36 
V g = 7.1 vT (m x v/p) 
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BU f bl ! p “ e J patlon de P th can be determined as a function of incoming fluid conditions 
using the following expressions for jet velocity decay (Reference 3-4) 


Vg/vj = 3. 09 Dj/AX 


(3-44) 


Do = 0.324 AX 


(3-45) 


where 


V 1 and D x - jet source velocity and diameter 


v 2 and Do = average jet velocity and diameter at a distance, AX, 
downstream of source 


Referring to Figure 3-16, jet conditions at the liquid surface relative 


to the source 


V L = 3. 09 v Dqo/H 


(3-46) 


D l = 0. 324 H 


(3-47) 


Also, jet conditions at X relative to re-entering conditions at the liquid surface are 
Vtv = 3. 09 Vt Dt /x 

L w l A (3-48) 

Substituting (3-46) and (3-47) into (3-43) gives 


v T .v = 3. 09 v D np /X 


(3-49) 


Solving (3-49) for X results in 


X 3. 09 v Dqj^/vj^ 


(3-50) 


Irom assumption 6, substituting 10 v^ = v . into (3-50) gives 


X — 0. 309 v V|^ 


(3-51) 


Finally, bubble stay-time can be determined by substituting (3-43) and (3-51) into 
t = X. vvj. which results in v 


t = 0. 102 Dq R v (dg) 
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Bubble stay-time Is plotted In Figure 3-13 as a function of momentum Inflow and 
acceleration environment. Note that an acceleration environment of 10-1 ,, lvlll 
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NOTES; 1. Vapor entrainment model from Figure 3-16 Is assumed. 

2. Dur *« stay-time, t-.l02 Dor v/(dg) 

3* Dubu.w diameter * . 41 cm (0. 16 inches), selected because it represent* 
maximum size (see Figure 3-15). 
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Figure 3-18. Fractional Volumetric Hold-Up of Vapor in Liquid Hydrogen 
as a Function of Liquid Momentum Inflow and Bubble Stay-Time for 
the Fluid Management System Tank 

result in bubble stay -times of five to nine seconds. Figure 3-19 provides plots of e 
versus momentum Inflow for a 50 percent liquid fill, and bubble stay-times of one, 
two and three seconds. These results should be used with caution, however, because 
the empirical data upon which Equation 3-35 is based was limited to Vg/VL? 20. Thus an 
uncertainty exists for ratios in excess of 20. Nevertheless, it is obvious that a mechanism 
does exist for the entrainment of large quantities of vapor during the tank fill process. The 
ensuing discussion will show that large quantities of entrained vapor are not needed to 
provide a near thermal quilibrium environment during fill. 

3. 3. 3. 5 Bubble Heat Transfer During Tank Fill. Equation 3-11 can be used to 
determine heat transfer from the dispersed bubbles to the liquid bulk during fill by 
Introducing the appropriate liquid properties. This equation becomes 


'P °l; 


2/3 r (m v 2 /V )Mr 1 1/4 

= 0. 163 — i 


(3-53) 


where h is now the heat transfer coefficient between bubble and liquid, and \'l is the 
liquid volume rather than tank volume. 
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Introducing liquid hydrogen properties at 275 kN/m 2 (40 psia) saturation, and solving 
for h, we have 


r • 2 \ l/ 4 

h = 812.1 Ft I > watt/m 2 - 0 K 


Now, the heat transfer rate for each bubble can be expressed as 


(3-54) 


Qb h V T B- T L> nd 2 , 

V B V B (tt/6) d 2 B L 


6h (T b -T l ) 


(3-55) 


where 


q_ = heat transfer rate from each bubble 
V B = bubble volume 


Tg = bubble vapor temperature 

T l = liquid bulk temperature 

Total heat transfer rate from the total dispersed vapor volume can be determined from 
Equation 3-55 by introducing the total number of vapor bubbles, n, 


Q t n Qg Qg 6 h _<VV. 

Combining Equations 3-54 and 3-56 gives 


(3-56) 


o 1/4 

— = 4372 (mv“/V L ) (Tg-T L )/d 


(3-57) 


Finally, 


Q = 4872 V (m v 2 /V l ) (T b - T L )/d 


(3-58> 


where 


total heat transfer rate from the entrained vapor 


Vg - total vapor volume entrained in liquid 
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Equation 3-57 Is plotted In Figure 3-20 for a range of ( and (inv 2 V^L Equation 3-58 
is plotted in Figure 3-21 as a function of ( and liquid percent fill of the fluid management 
system tank. Each figure assumes (T B -T L ) -.IK (0. 18R). Bubble diameter, d, is 
determined from Equation 3-33. 

3.3.3. t» Achieving Thermal Equilibrium During Tank Fill. An estimate can be made 
of how closely thermal equilibrium will be approached during tank fill by determining 
the ullage energy removal required to achieve equilibrium as a function of vapor to 
liquid temperature difference, and comparing to the heating rates of Figure 3-21. 

Ullage energy removal Is determined by employing the First Law of Thermodynamics 
for a control region that Includes the Initial ullage mass 
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where 


iEg ■ (Ujj nig)* - (Ug mg)j - change in ullage energy 
A El ' (Uj^ n^lo ~ (t?L nj lJi ~ change in liquid energy 
AQ - net heat transfer across the control region boundaries 


BUIUU.E STAY- TIME IN L*5UU>, SECONDS 
3.0 2.0 1.0 

... r — 1 p — ; r~T 

/ -/'\-A y"\ \ -i- -) 

/; /. : » y ; i j } 


Since the control region encompasses 

the' initial ullage mass only, m | 0 

ami m lV ■ nit m lV . Thus 
1 - -M 


• 3|' - V ■ > ^ -* 


/ i MOTS 

/ | t 


1, Tank VnUinttf o. .j 

If . I«1 ill] Ml 

/ i /'.. . J 3 ‘ • * v « 'V'V- v *''-'V‘’ v t.* i ‘ 

' / / !/ j j j »h«’rr ; 

! , j i 'll \ s K i* in.m M«ur».i-iT. .! The vapor ;uid liquid at condition 2 

I jr/y\ j j l''T!um*\iu^f.»r i are assumed to be in equilibrium 

■ •/ I ; i i 1111 . with the liquid bulk. Saturated 

| a,,... ; — ~ — | — j — » — j — j ■ j r hydrogen properties at 275. r> kN m- 

' h, , j ! . , | 1 , . i | | .-.i...... (10 psia) were selected. The product 

| 1 j ! ; : I ' ! j i M (lUgUglj was found to be approximately 

# 5 (30( constant with pressure and independent 

momentum inplow (riiv), iibm-«/»rc s ) °t temperature which meiuis that ,\Q 

is a function of Initial ullage pressure 
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Management System Tank and liquid Is directly proportional to 
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, \ r««k ri« and liquid, it Is, therefore, more conserva- 

, ...} tlve to assume that the vapor is Initially 
' \ saturated. The curves of Figure 3-22 are 

Y ’ solutions of equation 3-59 for \qi versus 

* ' • to.J : percent tank fill and (T|c I* i„) • where Ip 

^ [ \ , . , Y - j Is vapor temperature saturated at ullage 

\ j pressure. A comparison of Figures 3-31 
V ’ and 3-32 Indicates that equilibrium can be 

• achieved In less than one second. 
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Harure 3-22. Vapor Energy Removal Required to Achieve Thermal Equilibrium 
W Ith Liquid Hydrogen in Fluid Mangement System Tank 

considerable vapor entrainment will occur during fill, in spite of fact that the model 
is more applicable to a one-* environment than a very low-* environment, the belief is 
that considerable vapor entrainment in liquid will occur as a result of liquid inflow 
conditions Furthermore, the role of vapor buoyancy is such that bubble stav-tlmes 
will be substantially greater in ,-ero-g than In one-*. 

This leads to the conclusion that t hernial equilibrium is more likely to occur in zero-g 
than at normal gravity for the same propellant inflow conditions. 

Near-thermal equilibrium conditions should be experienced during low-* tank fill for an 
acceptable range of liquid inflow velocity and flowrate. The results of Figure .'1-21 „how 
that vapor to liquid heat transfer rates in excess of t?3. :* k\v pit) Wu sect cut be ~ 
expected to occur for inflow conditions of 0. 23 kg sec <0. 5 lb sec) and 2-1.4 m sec ,80 
ft, sec), which is equivalent to a fluid power input of CM. a Ug-m-’ sec”' ,:!200 lie -ft2 
sec' ). This ullage chilling condition is sufficient to achieve thermal equilibrium' in a 
traction ot a second. Although specific calculations have not been made for lower inflow 
conditions, a review of Figure 3-20 indicates that an incoming flowrate of 0. 0»a kg sec 
l . 1 lb, n sec) will still provide neap-thermal equilibrium fill conditions. 

:l. l VAPOR REMOVAL FROM .SCREEN DEVICE 

At the completion of propellant fill, vapor will be trapped within the screen device 
because procedures to prevent this are not known. Vapor can. however, be easflv 

removed trom the channels. A condensation process will be emploved for vanor 
removal. ' 



Vapor condensation begms when subcooling is Initiated, which will occur when the tank 
is pressurized with helium from 344.5 kN/m2 (50 psia) to 413. 4 kN,m2 m psia). 

There is no question that the trapped vapor will eventually condense. The onlv 
uncertainty is the time required to complete the process. 

3. 4. 1 BUBBLE CONDEN. TION MODEL. .An analysis was conducted to determine the 
time required for hydrogen bubbles to collapse. 

bubble condensation times, a bubble collapse program was written tor the 

anT ' pr0erun “ lcula, « ‘be collapse time of bubbles of different 
inttla 1 radii for various tank pressure differentials. The equations used In the program 

are based on the theoretical Investigation of Ftorschuetc and Chao, presented In ‘ 
Reference 3-a. The analysis Is applicable to any pressunmt as long as noncondensible 
gases do no. exist In the vapor bubble. The collapse rate Is controlled. prtaarUv b! 
eat transter rather than by liquid Inertia or surface tension. The analysis for heat 
transfer controlled bubble-collapse Is based on calculating ,h, conduction from the babble 
surface Into the liquid. The general equation Is 


1 -VT 


H 




where 


Y (nondimensional radius) = r/r 0 
T H (nondimensional time) = - j a 2 


•Ja - Jakob number, o C ^T/(p h fo .) 

L L V IS 

o thermal diffusivity of liquid 

t - time 
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°v 
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3 bubble radius 
' initial bubble radius at t - 0 
- liquid density 
vapor density 
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saturation temperature at final system pressure minus 
- liquid specific heat capacity 
heat of evaporation 


liquid temperature 






The numerical solution of the equation is presented in Figure 3-23, where nondimensional 
radius is plotted versus nondimensional time. 


Results of Equation 3-61 for hydrogen are presented in Figure 3-24. The curves plotted 
in Figure 3-24 are the times required for complete bubble collapse (Y = 0) for a given 

initial tank pressure, bubble rad? us. 
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Figure 3-23. Heat Transfer Controlled 
Bubble Collapse 


and tank pressure differential. It is 
noted that a bubble of 2. 54 cm (1. 0 
inch) initial radius will collapse within 
210 seconds. This dimension is 
substantially greater than that of the 
0. 64 cm (0.25 in) screen channel 
— width, which is believed to be the 

dimension for heat transfer controlled 
bubble collapse. The condensation of 

vapor contained within the screen 

channels will be completed within 
minutes after the system is pressurized. 
It is unlikely that this time to complete 
vapor condensation will be critical to 
system performance. Guidelines have 
not been established to show when 
pure liquid must be available within 
the screen channel for propellant 
transfer, nor has this study identified 
a need for early availability of pure 
liquid. 
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Figure 3-24. Collapse Time for Spherical Bubbles in Liquid Hydrogen 
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3.4„2 SUMMARY. Pressurization is an acceptable means of removing vapor from the 
screen channels, if the resultant condensation process is completed within approximately 
four to five minutes. 

3. 5 TANK REFILL 

Fluid management systems will require propellant refill capability in order to extend 
their useful life in space. The refill procedure must be capable of handling all operations 
from topping-off a tank with liquid, to a complete refill. Such capability is mandatory 
not only because propellant resupply can be a costly operation, but also because flexibility 
with orbital procedures must be maintained. 
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There are similarities between the initial fill of a fluid management system, and 
subsequent refills. Each requires pre-conditioning before liquid flow to the tank is 
Initiated. Once pre-conditioning is complete, propellant refill will be identical to initial 
propellant fill, which is discussed in detail in Section 3. 3. This section will deal only 
with the specific operations needed to pre-condition the fluid management system. 


Pre-conditioning will depend solely on the use of autogenous pressurization and helium 
pressurization. A helium pressurization system will require a pre-conditioning 
procedure because much of the helium must be vented prior to beginning liquid flow. 
The procedure may be complicated if helium enters the screen volume. Potential 
problems and solutions are discussed in Section 3.5. 2. 

Pre-conditioning will not be required for an autogenous pressurization system. Refill 
can be successfully implemented, as will be shown in Section 3. 5. 1, by flowing liquid 
hydrogen into the tank. 


3.5. 1 TANK REFILL (AUTOGENOUS). Refill of a system containing liquid and its own 
vapor as a pressurant is a straightforward operation. Liquid must be introduced at the 
correct vapor pressure, and fluid inflow power must be sufficiently high to assure near- 
thermodynamic equilibrium conditions. A determination has already been made that 
near-thermal equilibrium conditions will exist during fill. To be considered is the 
relationship between Initial and final tank fluid conditions, and entering liquid vapor 
pressure for a thermal equilibrium process. The inlet vapor pressure can be 
determined on the basis of the following First Law analysis Gi. a.-uV/U P\GE IS 

OF POOR OUAT.TTV 


dE L + dEg = hdm 

dE L = ( u L m L>2 “ < u L m L ) l = change 111 llquid ener Sy 
dE g = (Ugm g ) 2 - (Uging)! = change in vapor energy 


(3-62) 

(3-63) 

(3-64) 
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hdm = h (m T2 - m T1 ) = total energy of entering liquid 
m-j> = m L + m 

5 


(3-65) 


■" ,he *■* are 

assumptions with Equations 3-62 through 3-65 results if ' ^ = Ug ' C ° mbining these 


U L (m L2 * m Ll> f u g ( m g 2 " m gl ) = h (m T2 - m T1 ) 
Now, total mass within the tank can be expressed as 

m T = v t P L “ (or - 1) m g 


(3-66) 


(3-67) 


m T2 - m Tl = («-!) (m gl - m g2 ) 

where 

m T m L + m g = total mass of fluid in tank 
P L = liquid density 

P g = vapor density 
' * “L /P g 

Combining Equationd 3-66 through 3-68 we find that 

”L <m L2 • “LI* * U g< m g2- m hl> ' h <“ -D<m gl 



(3-68) 


(3-69) 


•Also, by adding and subtracting u L m gl and u L n, g , to the left side of Equation 3-69 
U L< m T2- m Tl> + u ev < m g2 “ m gi ) » h (or - 1) (m gl - m g2 ) (3 _ 7 , 

where 


u ev = (u g -u L ) = internal energy of evaporation 
Combining Equations 3-68 and 3-70 

u L (e - l)(m gl -m g2 ) -u Qv (m g2 -m gl ) = h fe - l)(m g2 -m j) 


(3-71) 
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Solving for the entering liquid enthalpy, h ORIGINAL PAGE IS 
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u ev 

\ a - 1 (3-72) 

Results are plotted for liquid hydrogen in Figure 3-25 as incoming liquid vapor pressure 
vorsus final liquid vapor pressure in the tank. It is interesting to note that this vapor 
pressure relationship is independent of liquid mass in the tank at the start of refill. 

The above results will simplify refill because the same entering liquid vapor pressure 
will be required regardless of initial liquid fill. For example, an entering liquid vapor 
pressure of 235. 6 kN/m2 (34. 2 psia) will be needed to maintain the tank propellant 
vapor pressure constant at 344. 5 kN/m 2 (50 psia). 

3.5.2 TANK REFILL (HELIUM PRESSURANT). Refill is complicated by the need 
to vent helium before propellants can be introduced. If helium is not vented, its 
partial pressure will become excessive as a liquid fill approaches 100 percent. 

This is illustrated in Figure 3-26, which shows the influence of initial liquid 
volume on final tank pressure, following refill to a 90 percent liquid volume. 

An initial helium partial pressure of 68. 9 kN/m 2 (10 psia) will increase tank 

pressure to 964. 6 kN/m 2 
(140 psia) as liquid is tanked 
from 10 percent fill to a 90 
percent fill condition. 

3. 5. 2.1 Tank Helium Vent. 
There are two methods for 
venting helium from the 
tank. The first is the 
transfer of liquid to a storage 
tank through the acquisition 
device until unusable resid- 
uals remain, and then venting 
the ullage until pressure has 
been reduced to a pre- 
determined level. The 
second method simply initiates 
ullage venting without first 
expelling liquid. The mass 
of vapor vented, to achieve 
a given helium pressure 
following refill, can be 
determined from Figure 
3-27. It can be seen that 
vent mass increases with 
a decrease liquid volume 


CONDHTONS 

1. TANK VOLUME! = 0. G2 M’* (22 FT^}. 

2. Lr-jctr va pop, ?resschs tn tan:; remains constant. 
‘H'Sr.HCDYN \:.!IC SQ'.TUBRIVM CONDITIONS EXIST AT 
LK? TEMPSFw^VSF. 

*• XC HELIUM PRESSURE NT. 

3 . da t.a is a ?r- l :c.a b l e t o any inttia l liquid ?:l l volume . 



( 30 ) (401 


LIQUID VAPOR PRESSURE IN TANK, kN/m 2 (pah) 

Figure 3-25. Entering Liquid Hydrogen Vapor 

Pressure Required to Maintain a Constant Liquid 
Vapor Pressure in Tank During Fill 
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C. NO HELIUM IS VENTED DI KING TANK FILL. 
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INlTiAl TANK FI LI , pcrciMi by volume 


CONDITIO NS 

1. tank volume - o. r 2 m 2 ^2 ft^l 

2. INITIAL LIQUID VAPOR PRESSURE - 3W.3 IcN/M 2 
(30 PSLU. 

3. IN IT LA PELTUM PARTIAL PRESSURE - *JS.D kNVM 2 
,;in psLAL 

t. FIN A L UET.IVM PARTIAL PRESSURE IS EASED U?n\ 
i I' 1 *’ V ’ L L . ' O E V' r ' L UK I E 70 L L i W IN ..I T N K E E FILL. 

; . V A PC R i o: L V IS "ENTER. 



Figure 3-26. Final Fluid Management Figure 3-27. Hydrogen Vent Mass Required During 
System Pressure for Propellant Refill System Blowdown to Provide Indicated Helium 

Without Prior Helium Vent Partial Pressure Following Propellant Refill 



at the start of venting. Figure 3-2S gives the tank pressure at the end of venting as a 
function of helium pressure following refill. This figure shows that the tank must be 
vented to lower pressure levels as the initial liquid volume is reduced. 

The data of Figures 3-27 and 3-23 appear to indicate that liquid should not be transferred 
before venting. These results, however,, ignore the possibility of two-phase venting in 
a low or zero-g environment. The likelihood of venting liquid can be reduced by 
decreasing vent flowrate. It is judged that an excessively long vent duration is required 
before venting can be conducted in a surface tension dominated environment. 

The preferred method for venting helium is to first transfer the propellant to a storage 
tank. Two-phase venting, for this condition, is minimized because any propellant mass 
lost overboard will be insignificant. 

The most convenient method for transferring propellant appears to be the use of the 
supply tank, which will service the system, as a fluid receiver. The following 
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CONDITIONS 

1. T.AXK VOLUME * 0. M* (22 FT'*). 

2. INITIAL LIQUID VAPOR PRESSURE ■ 344.5 kv/\f2 
(50 PS LA). 

3. DfITLAXi HELitT.I PAP.TIAL PRESSURE » «3.* k\V\'2 

1C ?« r \ ), “ ‘ ' k 

*• rec.AL HELIUM PRESSURE !S BASED U?OX *0 r 
ULLAGE VOLUME FOLLOWING TANK REFILL. 

3- VAPOR ONLY !S VENTED. 



Figure 3-28. M&nag oient System Pressure 
at End of Vent Required to Provide 
Indicated Helium Partial Pressure 
Following Propellant Refill 


procedure will be used: 

1. Chilldown the propellant transfer 
line connecting the system and 
supply propellant tank. 

2. Transfer all usable propellant 
from the system to the supply 
propellant tank. Liquid residual 
will be less than five percent by 
volume. 

The system will be ready for helium 
expulsion at the completion of Step 2. 
Helium will be expelled by performing 
a tank blowdown to reduced pressures. 

The above will not overly complicate 
a refill procedure and is preferable to 
introducing an extended vent duration 
coupled with the uncertainty of two- 
, 00 phase venting. 

Figure 3-28 indicates that system 
pre-conditioning will be satisfied 
once a tank blowdown Is performed 
to reduce pressure below 41. 3 kN/m 2 
(6'Opsia). According to the figure, 
refill to a 90 percent level will 

partial pressure 
Unfortunately, this 
entry to the screen 
requirement of providing 


, . , . maintain fi: 

below 68.9 kN/m2 (10 psia), which satisfies refill criteria. 

procedure does not address the potential problem of helium 
channel device, which would violate the primary system 
pure liquid flow to other systems. 


There is a possibility that liquid hydrogen evaporation during tank blowdown will cause 

helium T?TT * ^ AU »e taken to prevent 

lium entering because it cannot be expelled. Thus, the primary requirement, of 

provid r>g pure liquid flow to other systems, could be violated. The following 

describes the most likely mechanism for the flow of helium into the screen device: 

1* ° Uring tank blowdown to 41. 3 kN/m2 (6. 0 psia), the residual liquid will boil as 
its vapor pressure decays from 344. 5 kN/m2 (50 psiai to 34. 4 kN/m2 (5 psla) . 

2. The boiling process will result in considerable vapor generation « thin the screen 
volume, which will expel much of the contained liquid. 
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entry be prevented ^mott direct^thod^f ^ SCre6n d6ViCe demand that helium 
helium prior to initiating refill Two ancroa h Preventuig this is to ex P«l all available 
to a vacuum and, a tank pur». “ considered were, system evacuation 

consideration because of propellant freezing 3 , eiimlnated trom detailed 

below 6. 89 W/ m 2 (1 . o J t Tn* £Z °^orZt y TOg ' n Wffl ,re ~ e - P-sures 
cannot be ignored as tank pressure fs reduced J 1 ! ZUlg W “ hfa the ve ”‘ s y stem 

purge procedure was adopted. a vacmun * Consequently, a tank 


concentrations wn^chwm nou^rfweT, 1“™ 0a “ be dIluted t0 

repeated charge and vent cycles wifh th een device ' s functioning by 

than 0.25 cm fo. 1 hiXgZ * heH ““ ™ 

that should not interfere with the effective fun Z"" 3 represeMs a total volume 
to a total helium mass of 0. oo« k^ o ooU™ Ttu * ^ ^ ™ S S ^»alen. 
expended to achieve this total helirnn mass wfs h ^ ““T* of gaseous hydrogen 

maSS ’ Was based the following procedure: 

1 an. 


X t 7 fT ae tok * 1U * .o 

Ihm) c( helium remains in the tank. ‘ Lnde ‘ the3e conditions 0. 0 03 kg (0 . 0183 


2 . 
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Step 2 w„l he repeated until no more than 0. 0013 kg ,0. 0020 Uy of helImn rcmaUls . 

OM/h O ftfS 1 
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Figure 3-29 shows the results of char°ine tossnv/ 2„ n 
(20 psia). Helium mass is reduced to an acceotjihl* 1 PSla) 31,(3 t0 137 * 9 

cycles, and a consumption of 0. 25 kg (0 56 lh \ a n Tn V/ith tW ° charge and vent 
This figure indicates that sufficient helium dlfei d °* u° Rg ( °’ 66 lbm) res P e ctively. 
period of time, with only a smaiuoss of v °“ be a ° hleved a relatively short 

tank to 68.9 kN/m2 ( io pS ia) is selected b^aus^thehTlf CyCl * that charges the 

effective. 0 Decause the helium dilution process is more 


ZpZoToZTan^^^^^^ will he initiated np. 

relatively high vapor pressure due ‘Z- “ * 

hydrogen temperature. Figure 3-11 shows that th! 1 ? J , evacuatad tank at liquid 
exceed 289 kN/m^ (42 psia) if a final vapor pres^f^ 
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Figure 3-29. 


Fluid Management System Helium Residual Following GH 2 Purge 


e^e^Mced^ur^g^filHor ^°acc'epiabf/' t ^ CTma f\ e ^ U ^ i ^ rilm will be 

flowrate. The discussion of Section 3 3 f hqmd inflow velocity and mass 

low as 0. 045 kg/sec (0. 1 lb™ 'sec) will nr a* ndlCates that enterin S flowrates as 
refill. ( bm/ S6C) ulU provide near-equilibrium conditions during 

c^pUonTpr^peU^t teflli° vapor wmb D ?’ ICE F ° LL0WING REF ^L. At the 
process of vapor removal following propellant *°T ^ 

and the selected procedure is eoualW nn „ ,, k , U discussed in Section 3.4, 
procedure will be to pressurise the Lo™iwf *? COnditions follow ing refill. The 
helium, which will subcool the liquid bv n^sldf' 't 2 ^ 

Section 3. i that trapped vapor will be ennrW aT 10 P id> * 1 Was indicated in 
sam« time 1, applicable to the poet-refill eohdeneetta^rS,"^ SeC ° ndS - ™ S 
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a helium pressurization system, however, because much of the helium must be vented 
before refill is begun. The selected procedure for this case will be: 


1* Transfer all usable propellant to the . -supply tank. 
2. Vent the tank down to 6. 89 kN/m 2 (1. 0 psia). 


3. Conduct two cycles where the tank is charged with hydrogen vapor to 68. 9 kN/m 2 1 • 

(10 psia) and vented down to 6. 39 kN/m 2 (1. 0 psia). A total of 0. 25 kg (0. 56 I b. ) j : ' ' 

hydrogen vapor will be consumed in the process of diluting f he helium to an * i i I 

acceptably low level. ;> T. ) 


The tank is now ready for liquid refill. j ! : 

/ j ' 

Liquid refill for the system with autogenous pressurization will be conducted with liquid I 

entering at a vapor pressure of 235. 6 kN/m 2 (34. 2 psia). This inflow condition will ' 

provide a final vapor pressure of 344. 5 kN/m 2 (50 psia). It was estimated in Section I • ■ 

3. 3, that near-thermal equilibrium condition will be maintained during initial fill, > ; • 

at liquid flowrates as low as 0. 0<t5 kg/sec (0. 1 U^/sec). Thermal equilibrium for I 

liquid refill is expected at the same flowrate condition 1 

I i | 

Liquid refill conditions for the system with helium pressurization will be identical to 1 i ’’ 

those described above, except that entering liquid vapor pressure will be approximately ! i ' 

289 kN/m 2 (42 psia). Vapor will be trapped within the screen device at the completion ; j 

of refill. This vapor will be condensed within approximately 210 seconds after the ; j ' - 

system is pressurized to its operational level of 413.4 kN/m 2 (60 psia). j j j t 

j i 

I , I 
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FILLING ANALYSIS FOR SELECTED CONCEPT (EARTH STORABLES) 


The same fluid management system described in Section 3 was evaluated for filling with 
N 2 0 4 under orbital conditions. This system configuration reflects design requirements 
imposed by the problems of liquid hydrogen storage in orbit for long durations. There 
are design problems peculiar to the storage of N 2 0 4 But these problems , in -eneral 

contract C ° mpr ' tibility “ d ’ are beyond the scope of tWs 


System fill with N 2 0 4 in an orbital environment poses fewer problems than liquid 
hydrogen because of differences between an earth storable and crvogen system- 
temperature environment and operating pressure level. System temperatures will 
not van- significantly from N’ 20 4 temperatures at any time during a fill or refill 
operation. Consequently, neither prechill nor chill will be required to precede tank 


The management system operating pressure of 2067 kN/m? (300 psia) is considerablv 
higher than for liquid hydrogen. This higher pressure level will simplify tank fill ' 

™™f, S WU1 be maintainea belc '" «»• normal operating levels without 
difficulty. This is illustrated in Figure 4-1 which gives the maximum tank presr.me 

as a function of percent liquid fill. Isentropic compression of the ullage is assumed 
which meaiis that heat exchange with the liquid and tank walls is assumed to be rero ’ 
Note that 90 percent fill can be effected without exceeding 1378 kX.-m- ( 200 psia) 
pressure. Inlike cryogenics this extremely conservative approach can be used because 
results are acceptably low. 


One problem that cryogenics and earth storables have in common is that an undefined 
(or poorly deiined) liquid-vapor distribution will exist during low-gravity tank fill 
This inability to determine vapor location will make it difficult both to assure that' io 
liquid is vented during tank fill, and to assure that vapor is not trapped with:-, me 
screen channel device. 


A second problem common to an earth storable and cryogenic management system «s 
the need to vent helium before propellant refill can be initiated. Helium venting must 
be conducted with care under orbital conditions because liquid mav also be vented )n 

addition, helium may enter the screen device during this period, unless precautions 
are taken. 
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CONDITIONS 

I- TANK VOLUME 0. (>2 M-> (22 FT 3 ). 
INITIAL VAPOR TEM Tj . AND PRESSURE 
“ 30 °K (540R) AND Uil kN/M- ( 19 PSIA). 
NO HELIUM. 

3. THE VAPOR IS COMPRESSED 3 

ISENTROPICA LL Y DURING FILL. I 
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PERCENT LIQUID FILL 


Figure 4-1. Maximum Management System Pressure During Fill 
With N 2 0.{ 

4.1 THERMAL EQUILIBRIUM TANK FILL 

The single requirement for tank fill is that acceptably low pressures be maintained 
during the process. Thermal equ, librium fill represents the minimum propellant tank 
pressure condition and, therefore, should be maintained during fill. The thermal 
equilibrium tank till relationship derived for liquid hydrogen < Equation 3-2S) applies 
.is we 1 to This equation is given below, 

h L u g 2 m R.> n 'l.o * U L „ • .1 u,.. m in mi ,1-1) 





a enthalpy of entering liquid final vapor 
U £2 m ?2 * internal energy 

m L, = final liquid mass 

V = management system energy change during tank fill 
u L .2 35 final liquid internal energy 


thir^f I.KW ^ “ L2 U S= to ll,,uld vapor P>-« S ur. and temperature under 

LToZ ° 0n ? m ” nS \ Thus u L 2 “« »g2 «e known once final vapor pressure 

la specified. Au w can also be determined, when Initial and final temperatures are 

8 yen. Finally. h L (which is a function of entering liquid vapor pressure) can be 
determined for a desired liquid fill condition. 


Equation 4-1 is summarized for N 2 0 4 In Figure 4-2 which gives entering liquid vapor 

n7« IT* “.f funCtlon of tol,lal ta " k temperature and final tanked liquid vapor pressure. 
I. is in ereattag to note lhat the effect of Initial wall temperature upon final tank pressure 
is negligible. This data illustrates that the N 2 0 4 thermal mass will so dominate tank 


CONDITION'S 

1. TANK VOLUME = 0.62 as 3 (22 ft3) 

*• tank mass = 13. si kg ^30 lb^j, 2219 aluminum 

3. TANK PRESSURE IS 0. 0 kN/nr t 0. 0 PSlA) PRIOR TO FELL 

NO VENTING DURING FILL ORIGINAL PAGE IS 

OF POOR QUALITY 
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4.2 REFILL ANALYSIS 

There are two requirements that must be satisfied during system refill On. , - ,u , 
he Hum must not eater the screen channel device at any tfm. Tte «cL 
is that propellant tank Dressur* k i 5 The becond requirement 

(350 psial P ma '" be ‘° W the ““'“um allowable of 2112 KN/m2 

It is assumed that tank pressure will onr 7 i KT o 

z *. — ^r g 
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to evaporate liquid or sublime frozen propellant. °° ntents " ith warm va P° r In ord ^ 
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1. The partial pressure of helium and N 2 O 4 vapor will decrease as tank pressure is 
decreased. N 2 O 4 evaporation will begin at the screen surface once liquid vapor 
pressure exceeds the partial pressure of N 2 O 4 vapor. 



2. Liquid will be drawn into the capillary device to replace the liquid lost due to J ' ] 

surface evaporation at the screen. | ; 

• i 

i ! 

3. The liquid bulk cannot boil while pressure exceeds liquid vapor pressure. j 

Consequently, tank vent will be terminated at a pressure greater than 131 kN/m2 | 

(19 psia) to assure that the capillary device remains filled. I , j 


Note: The only evaporation (or boiling) that can occur is at a liquid-vapor 
interface, and only if vapor pressure at the surface exceeds partial pressure 
above the surface. 


A determination has been made of the amount of venting required prior to tank refill. 
Figure 4-3 shows the final helium pressure following refill as a function of tank 
pressure at the end of vent, and the percent liquid in the tank prior to venting. It 
is noted that a final helium pressure less than 1936 kN/m^ (281 psia) will be realized 
without initially venting the tank below the N 2 O 4 boiling pressure. Figure 4-4 
indicates that approximately 3. 2 kg (7 lbm) to 4. 5 kg (10 lbjjj) of N 2 O 4 vapor will be 
vented during this process. 
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CONDITIONS 

1. TANK VOLUME = 0. 02 m3 ( 22 ft 3 , 

2. INITIAL N. : 0 4 VAPOR PRESSURE 
- 131 kX/ni2 (19 psiai 

3. INITIAL HELIUM PARTIAL 
PRESSURE = 193i> kN/m- (2S1 
psia) 

». INITIAL FLUID TEMPERATURE 
= 30 OK (540R) 

5. FINAL HELIUM PRESSURE IS 
BASED UPON 10~ ULLAGE 
VOLUME FOLLOWING TANK 
REFILL 

G. NO LIQUID IS VENTED. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


(10) (03) (30) (40) (50) (60) 

TANK PRESSURE AT END OF TANK VENT, ’icN/m2 (psia) 

Figure 4-3. Management System Pressure at End of Vent Required to Provide Indicated 
Helium Pressure Following N 0 O 4 Refill 
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quantities are present. A requirement exists hmv ^ k ® y t0 ° CCUr when 3m aller 
between the liquid bulk and the screen channel’ liquid' ^Thi^ maintain commu ^ations 
to assure that the bulk replaces losses due to COI ™'*ication is needed 

Sufficient liquid must be present in th* m i,? rf evaporation at the screen. 

mairtained. it was judged that a ten percent liJSd w?*d COmmunications W «1 be 
represents an acceptable compromise to this conflict " ^ Stait ° f helium vent 


4.2. l.i Propellants Consumed During Helium vw o * 

at the loss of some propellant and helium pressurant’ Fi^T 1 7? 1 be accom P lis bed 
matelv 3. 3 kg (3. 3 Ibm) of N 0 O 4 vapor will be e\Delli>H h 4 ~ 4 Sh ° WS that a PP«>xi- 
1 * 7 k? (3 * 7 lb m) of helium will also be vented The^ #1 ^ ^ ea ° h refiU P rocedure : 
regardless of the selected refill procedure Thp f ?, ?u quantities "<U be lost 
vented before refill can be initiated and / • ? ^ S thnt the helium must be 
be vented in the process. ’ mLX ‘ ure of heUu ™ and propellant vapor will 
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4 ' 3 VAP ° E BEM0VAI ' SCREEN CHANNEL DEVICE 


Upon completion of prope llan t fill n-r r o«n 

device, because procedure- to prevent this are 00^ ** trapped wtthln the screen 
removed from tbe channel CZlZiT tITs' Vapor «*• *»~». be 

,u be emp,oyed f ° r -««- - *» eCnie b ;r s ;x;. s 5 r prooess 

level of 2067 wlS POO^S? 1 "^ Tr^I? SySt6m *" pressuriz€d to its operating 
immersed, or surrounded by, subcooled TrorJlTT WlU condense due to being 
t e heat transfer controlled bubble collapse D ro ^ ^ general ^“ation describing 

equation is also applicable to the ccTden Z7on 7Zo That* 

3-61 for N2O4 properties are given in Fi™7 /I 2 4 Vap ° r - Results of Equation 
considerably shorter than for UquTd t d rT tl Bubble C ° lla P se «“es are 

sec 4 T (1 ; ° iD) taitial radius * An N2O4 vapor bubble T *? °° mparin e bubb1 ^ of 
seconds, when immersed in liquid subcnnl^fh 7* 11 colla P se in less than 10 

Mrogea bubble will 0011^.^1?^? ? kN/ “ 2 < 4 ° P*d>i "hereas a 

As discussed in Section 3. 5 the timp tn n 6r ltS system operating pressure 
device will not exceed the ttoe for^ 2 Within the screen ’ 

a 2. 54 cm (1. 0 in) radius bubble to collapse. 



Fluid management system refill will be iri r 1 ^ ^°° R QUALITY 

has been completed. For the case of initial" tai fm ^ flU ° nCe Pre-condition ing 
required. Liquid will simply be introdueeH , , i ’ n ° pre " con ditioning will be 
-d velocity to maintain a high flowrate 

analyses were not conducted for N 2 CL it is mM ®°“ d ? tona - Although specific 
0^15 kg/sec (0 . 7 lbm/sec) wjU ^4“ ^ e S «”“ »W- flowrates as low as 
process. Venting will not be required during tank fill. P ressu res during the fill 
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Although no precautions will be taken to keep vapor from entering the screen device, 
the collapse time is so rapid for N 2 O 4 bubbles, that the device is expected to be filled 
with liquid at the completion of tank fill. This condition will exist because fill will 
require several minutes to complete, during which time the propellant will be subcooied 
by more than 276 kN/m 2 (40 psi). As previously indicated, vapor collapse times for 
this subcooled condition will be substantially lower than the time to complete tank fill. 

Pre-conditioning is required for system refill because much of the helium must be 
vented before refill is begun. The selected procedure for this case will be: 


1. Transfer liquid to or from the tank in order to provide a ten percent fill condition. 



2. Vent the tank down to 206.7 kN/m 2 (30 psia). 

The screen channel device will remain filled with liquid during this time, as all 
evaporation at the screen will be replaced by the tank liquid. A total of 3 . 8 kg ( 8 . 3 lbm) 
of N 2 O 4 and 1. 7 kg (3. 7 lbm) of helium will be vented during each refill. These fluid 
quantities will be vented regardless of the selected vent procedure. 

Propellant flow conditions for tank refill will be the same as those identified above 
for the initial fill condition. Tank pressures will be maintained acceptably low during 
this process so that venting will not be required. Furthermore, the screen device will 
remain filled with liquid throughout this time period. 



MODELLING OF PROPELLANT FILL PROCESSES 


Fluid management system fill and refill is examined in this section to determine what 
experimentation is required to demonstrate this capability in a low g environment. 
Emphasis is placed on identifying normal gravity tests that can be conducted to model 
low gravity filling behavior. This approach was taken because of the ease with which 
ground based tests can be conducted compared to the more complex, and considerably 
more expensive, space based tests. The modelling and scaling techniques employed, 
however, are also applicable to space based tests. 

Based on the analysis performed in Sections a. 4 and 5 the important phases in filling 
cryogenic and noncryogenic tanks containing "total control" type acquisition devices have 
been Identified. Refilling cryogenic tanks containing capillary devices that are partially 
full demonstrates the complete range of possible phenomena that must be handled. 
Removal of the liquid from the tank must be accomplished prior to the removal of helium 
from the tank. As discussed, removing the helium by purging and venting will be one of 
the most difficult phases of the refilling process. Once the helium is effectively removed 
refilling is similar to the initial filling of the tank. 

Prechill of the tank is accomplished by inflowing liquid or cold vapor to cool the tank to 
the desired level. When the fluid in the tank has reached predetermined pressure and 
temperature conditions the tank is vented. Venting only occurs when vapor exists in 
the tank. The prechill/vent sequence cools the tank sufficiently to allow the tank chill 
process to subsequently be carried out with a locked up tank. This eliminates the need 
to control the liquid/vapor Interface for venting purposes. 

After tank chill has occurred, the filling process is identical for cryogenic and non- 
cryogenic tanks and efficient filling will be accomplished by maintaining thermodynamic 
equilibrium. Spraying liquid into the tank at a high flowrate and velocity, provides the 
fluid agitation to accomplish a low pressure rise condition, and maximize the amount 
of liquid that can be put in the tank without exceeding a given tank pressure. 

At this point the tank is full but the capillary device is likely to contain trapped vapor. 
This vapor will be condensible, if the helium removal has been successful. The vapor 
will be condensed by subcooling the trapped vapor with pressurant. In cases where 
warm vapor is used for pressurization, the pressurization should occur just prior to the 
end-use propellant outflow from the tank in order to minimize any mixing between the 
liquid and vapor that could cause tank pressure reduction and subsequent bulk boiling in 
the capillary device. 




5. 1 PROCESS SCREENING FOR MODELLING 

A thermodynamic model was developed in Sections 3 and 4 for each of the processes 
identified above. These analytical models were combined to provide a complete 
description of system fill and refill, and serve as the primary basis for identifying 
problem areas and developing procedures for low-g fill. In several instances th. 
models provide an accurate representation of the physical process. Experimentation 
will not be required for these processes. Additional data will be required for the 
other processes, however, before sufficient confidence can be placed in procedures 
for low-g fill and refill. 

The first task, therefore, is to conduct a screening evaluation for placing processes 
in the following categories: 

1. Process is subject to rigorous analysis. 

2. Process is adequately defined, but uncertainties exist regarding analytical models 
or the ability to extrapolate to low-g. 

3. Process is not sufficiently well known and, therefore, extrapolation to a low-g 
environment Is unacceptable. 

Processes that belong in the first category will not be considered for experimentation 
because the analytical models are sufficiently accurate to enable evaluation. Similarly, 
the third category of processes will not be subject to experimentation because in- 
sufficient information is available to identify an acceptable test plan. Only processes 
in the second categoryjidll-be^valuated to determine modelling and scaling eauations 
from which norjnaTgravity tests can be conducted and results extrapolated to a low-g 
environment; 


PROCESSES SUBJECT TO RIGOROUS ANALYSIS. The only processes in this 
category are tank vent and purge for a system which must expel helium. 

5. 1. 1. 1 Purge of Cryogenic System Prior to Refill. The ullage vapor of a nearly 
empty propellant tank will contain hydrogen vapor at 344. 5 kN/m 2 (50 psla) and helium 
at 68. 9 kN/m 2 (10 psia). The tank will be vented down to 6. 89 kN/m 2 (1. 0 psia) 
preparatory to a purge that will dilute helium concentration to an acceptable level. 

Two charge and vent cycles will be conducted using warm hydrogen vapor; the charge 
period will increase tank pressure to 63. 9 kN/m 2 (10 psia) and the vent period will 
reduce pressure to 6. 89 kN/m 2 (1. 0 psia). 

Analyses performed to describe the above are sufficiently accurate to obviate the need 
for experimentation. There is nothing in the charge and vent phenomena that is 
influenced by gravity environment. It is likely that any purge tests conducted will be 
useful in establishing a more effective procedure; but, such data is not considered to 
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be technology oriented and is, therefore, not applicable to a test plan, 

5. 1. 1. 2 Earth Storable System Vent Preceding Refill. It will be necessary to expel 
much of the helium contained within the tank before N 0 O 4 refill can be Implemented. 

The analysis of Section 4. 2 indicated that sufficient helium will be expelled in the 
process of venting the tank from its 'nitial pressure of 2067 kN/m 2 (300 psia) to 206. 7 
kN/m“ (30 psia). Bulk boiling of N 0 O 4 will not occur during the vent process because 
tank pressure will remain above the N 0 O 4 vapor pressure of 131 kN/m 2 (19 psia). The 
small liquid quantity contained within the tank will remain relatively free of disturbances 
as long as tank vent is conducted slowly. 

5. 1.2 PROCESSES REQUIRING EXPERIMENTATION. The three processes requiring 
experimentation are tank prechill, tank fill, and vapor bubble condensation. Each was 
analytically modelled (see Sections 3 and 4) in an effort to Identify acceptable system 
fill and refill procedures. The models do require verification, however, because the 
equations used represent extrapolations of existing empirical data. These processes 
will be discussed in greater detail in Section 5. 2 . 

5. 1.3 PROCESSES NOT ADEQUATELY DEFINED. One process that may be in this 
category is tank chill. Boiling will occur during tank chill as liquid spray Impacts the 
hot tank walls. This boiling phenomenon is rather complex due to the transient naturt 
of a vapor film being created and destroyed with each droplet impact. Equation 3-13 
describes the phenomenon of liquid splattering on a hot surface, and indicates that 
heat transfer will be independent of acceleration environment. However, the equation 
is considered to be unsubstantiated because there is no variable acceleration data 
available to confirm this independence. 

Rather than grapple with the boiling phenomenon uncertainty, the preferred alternative 
was to identify a procedural change to circumvent any uncertainty. In this case the 
decision was made to establish a prechill process to eliminate any concern with 
excessive pressures occurring due to boiling during tank chill. This was done by 
allowing prechill to remove sufficient tank energy that peak pressure would be 
acceptably low during tank chill. 

5.2 MODELLING OF SELECTED PROCESSES 

Modelling or scaling analyses will be used to determine the feasibility, and subsequently 
the conditions and configurations, of ground based model testing for providing data 
useful to determining the filling characteristics of an earth storable and cryogenic 
liquid fluid management system. 

5.2 . 1 TANK PRECHILL. The initial management system fill procedure in low-g will 
require a prechill process to reduce tank tempe<ature to an acceptably low level. 

Tank prechill will be accomplished by introducing low temperature hydrogen vapor at 
a sufficiently high velocity that forced convection heat transfer will dominate the heat 
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where 


h = heat transfer coefficient 

(mv 2 /V) = fluid power input per unit volume 
P» Cp, NpR and u are fluid properties 
n, m, A are empirical coefficients where 


n = 2/3 
m = 1/4 
A — 0. 163 
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Equation 5-1 will be In the magnitude of the empirical coefficients. The coefficient most 
likely to change is A; coefficients n and m are not expected to change. 

Ground-based tests for determining the validity of Equation 5-1 will be applicable only 
a free convection heat transfer is dominated by a forced convection process. It at... 
that the forced convection - 

the fluid input cower term, mv2. Thus, a series 01 precmu 
conducted, for a range of mv- that will provide forcrt conveotton ^mmmc ^ 
feasibility of this approach has been demonstrated with the aid ot the Appendix B 
subroutine where prechill of the Centaur-OSS model test tank was modelled Analyse 
results are plotted in Figure 5-1. It is noted that forced 

about three times greater than free convection at entering fluid conditions ot 0. 0009 
k"“orC/sec, 152 ra /sec ,500 ft/sec, . This inflow condition ,s feastbie 

for ground test purposes. 

nail Test Parameters. Tests will be conducted for a range of fluid power input 
values mv- using gaseous hydrogen and nitrogen flow into the Centaur-OOS tes an '. 
Hea^ transfer rates l the vapor from the hot tank will be determined using gas cmpera- 
ture and pressure, and tank temperature. The resultant heat transfer coefficien 
will be calculated. Empirical coefficients will be selected which provides the e* 
match for Equation 5-1 with test data. 

500 TANK FILL. The Intent of the tank fill process will be to create conditions 
conducive to attaining near-thermal equilibrium. These conditions will < . achteve 
bv introducing liquid into the tank through spray nozzles, which resulting spray 
provide the hfgb heat transfer rates needed to attain near-thermal equtlibnum 
conditions. As tank fill continues, the internal tank fluid environment will ctange 
from one where liquid droplets reside within the ullage volume to one where lapor 
todies are entrained within a liquid bulk. The transition from heat transter dominated 
hv llouid droplets to heat transfer dominated by vapor bubbles is expected to occur in 
2.' X of C «« «o 50 percent liquid fill. As propeilani tank «U continues he 

:ir - 

Consequently, this mechanism will be analyzed to determine what modelling 
scaling techniques will be applied toward the development ot a te^t plan. 


I 



Figure 5-1. Predicted Forced Con\ection and Free Convection Heat 
Transfer Coefficients During Normal Gravity Test Tank 
Prechill Test 


The applicable expressions for determining heat transfer between the dispersed vapor 
phase and its surrounding liquid were developed in Section 3. These equations are 

given below. 

q t = 4872 Vg (mv 2 /V L ) l/4 <T B ' T L )/d (from Equation 3-58) (5-2) 


0.6 0.5 

d _ h 134 -^ £ + 0. 09 (from Equation 3-33) 

y 0 " 4 (mv 2 /V L )°’ 4 

(5-3) 

v g- V 

(5-4) 

(mv/p L ) (from Equation 3-38) 

(5-5) 

t ^ v (D or /'dg) (from Equation 3-52) 

(5-6) 

« • V (V g + V L } 

(5-7) 


Combining Equations 5-4 through 5-6 
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All variables influencing the thermodynamics of propell^t .onk fill are contained in 
Equations 4-2, 5-3, 5-7 and 5-S. The task now is to establish a relationship. An 
inspection of these equations indicates that rigorous modelling techniques cannot be 
employed to establish equivalence between normal gravity and lov-g fill. This is due 
to the unknown influence of accelerat.m environment upon vapor entrainment, and 
because vapor entrainment is an important variable in the determination of Q T and d. 
The following discussion will show, however, that normal gravity tests will provide 
data useful for evaluates and verifying low-g fill procedures. 

The first step in this modelling assessment is to modify the above equations to reflect 
the ratio of ' onditions between model tests and orbital application. Let 

Q t * = Qtm/Q tp . d* = djn/dp, V g * = Vgm/Vgp» 


where 


m = one-g model test condition 


p = low-g prototype test condition 


Thus, the equations become 


Q t * = V g * (m v“/V) + 1/4 (T r - T L f / d* 

(5-9) 

d* = f* /(mv-/V) 4 

(5-10) 

p> + 

V g * = mv- (D or */d* g*) 

(5-11) 

c* = V g * / <V g +V L > 

(5-12) 


Fluid properties do not appear in these equations because it is assumed that the same 
fluid is employed for both the model and prototype test conditions. 

Equation 5-11 is the key to this assessment because its Influence is present in the 
remaining equations. More important, it is subject to the greatest uncertainty because 
of the acceleration influence, g*. This equation indicates th.T Vg* will decrease as g* 
increases, furthermore, g* -> > 1 when a one-g °nvirorment is compared to an orbital 
environment. Because the influence of the other variables of Equation 5-11 will 
obviously bo small compared to g*. V g - < 1.0. If this is so, then it follows that. 



•\ physical description of the above is given by describing changes to the internal tank 
fluid behavior tor a test where all parameters are kept constant with the exception of 
acceleration environment, which is gradually being decreased. 

1. A reduction in acceleration environment will increase the stay-time of entrained 
vapor due to a reduction in bubble rise velocity. 

-• This increased stay -time will result in a greater steady state volume of entrained 
vapor, because the rate at which vapor is entrained into the* liquid does not change 

The greater volume of entrained vapor will increase i *> d-v:ipor surface. 

1* Finally. ino combination *>t increased heat transfer area and time will result in 
an increased heat transfer rate between liquid and entrained vapor. 

It is •oneluded that ;he influence of a lo*\-g eavt ronnient will i** to increase the 
likel ikood that thermal equilibrium conditions . *1 e\i ;t during the lank till process. 

!t can, t heretore, be uateii that till tests conducted in normal gravitv represent a 
more sever-* environmem fot maintaining acceptably low tank pressures. 

It will still 1 h' extremely useful to conduct normal gravity tests if the results are that 
near-t-quibbrium cond.itiens will be m attained during !'i 1 1 . Such test results will be 
interpreted as reflecting a greater deviation from thermal equilibrium than will occur 
in a iow-g environment. With this limitation op the interpretation of test results, a 
Haw test parameter can be sheeted. Equations >-:> and .*-10 shew that the only other 
parameter to influence heat transfer of the dispersed vapor phase during tank fill is 
mv~ Y, the fluid power input per uni! volume, '’his fluid power input parameter 
attects seat iraustei eoettieient as well as vapor bubble diameter. It is recommended 
that irir ' be employed for scaling ground tent results to a lovv-g environment. 

Figure illustrates expected results fot* a series of tank fill tests. It is reasonable 
assume that tank pressure will approach equilibrium pressure as the fluid power 
input parameter is Increased. As indicated bv Figure >-.! a critical parameter value 
exists for which equilibrium will be maintained during tank fill. Anv additional 
increase in tin* parameter value will not further decrease propellant tank pressure. 

I'he test objective will be to obtain a relationship between the fluid power input 
parameter and* tank pressure, which is a measure of fluid equilibrium conditions, of 
special Interest w ill be (dent ification of the minimum parameter value for which 
equilibrium will occur. IVst.* will lie conducted using liquid hydrogen for the* cr\ wen 
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Figure 5-2. Expects Trend ot Tank Pressure History 
During Fill Process 


system tests, and Freon MF tor the earth storable system tests. 


rstem * * 

5.2.3 VAPOR CONDENSATION. At^the it system has a 

provlde^vapor-free liquid to other systems. Consequent vapor must 
be removed from within the screen device. 


A condensation process will * 

Tl : Uas based upon conduction heat transter controlled 

bubble collapse and Is expressed as 


H 


(from Equation 3-»>H 
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5 _X 3 and unsteady state heat conduction to or from 


There is a parallel between Equ^.on^ 
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Figure 5-4 Illustrates the process of vapor collapse due to condensation. It is assumed 
that conduction of heat from the vapor Interface and from the interface to the liquid bulk 
will control the rate of condensation. Liquid will enter the screened channels as the 
vapor volume is reduced by condensation. Although flow across the screen will be 
uniformly distributed, the c hann el shape dictates that the liquid will collect at ends of 
the channel. Thus the liquid-gas interfacial surface area will be reduced as vapor 
collapse progresses while the vapor width remains constant. 

There is no available analytical model for vapor collapse of the indicated condition. 
However, it is reasonable to assume that the Fourier Modulus, c*t/L-, is the key 
parameter in any evaluation, where 2L is the channel width. Because the model 
indicates that vapor width will remain constant, it is likely that the relationship of 
energy removed to dimensionless time will be like that for a slab. This relationship 
is given in Figure 5-4. 

This model suggests that ground tests must be performed in such a manner that free 
convection effects are less significant than conduction. Since free convection effects 
cannot be eliminated from the propellant fill tests, a screened device will not be 
included as part of the test apparatus. Instead, a series of tests will be conducted with 
liquid hydrogen and Freon MF which should be adequate to verify an analytical model 
both for earth-storable and cryogenic propellants. An evaluation of test results will show 
whether the model will be more similar to unsteady state conduction heat transfer in a 
slab, or to a collapsing volume model. 
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TEST PLAN - FILLING OF ORBITAL FLUID MANAGEMENT SYSTEMS 


The discussion of Section 5 identifies three processes of fluid management system fill and 
refill, which require experimentation to verify the analytical techniques developed in 
Sections 3 and 4. These processes are; tank prechill, tank fill, and condensation of trapped 
vapor. A modelling analysis has identified the dominant heat transfer and fluid mechanic 
parameter which controls the thermodynamic '.ate processes during prechill and fill. This 
parameter is the fluid power input' per unit volume or, mv 2 /V, and it is a measure of the 
intensity of fluid agitation within the pvopeilant tank. A determination was made that the 
vapor condensation process will be controlled by conduction heat transfer. 

The parameter, mv / V, will be used to determine forced convection heat transfer between 
the entering propellant vapor and hot tank walls during prechill. The basic form of the 
applicable heat transfer coefficient expression is given by Equation 3-11. Its empirical 
coefficients are uncertain because the expression was developed for a somewhat different 
set of conditions than will be experienced during prechill. Consequently, experiments 
will be conducted to identify a set of coefficients applicable to the prechill process. 

The parameter, mv 2 /V, is the dominant variable in identifying the internal tank fluid 
condition during tank fill/refill. It determines the volume fraction of vapor entrained 
in liquid, and the bubble size population. Thus, mv 2 /V is the key to determining heat 
transfer rates between entrained vapor and liquid during tank fill. Experiments will be 
conducted to Identify the threshold value of mv2/v for which thermal equilibrium will occur. 

The condensation rate of vapor trapped within the screen device will be a function of 
vapor volume and shape. It is obvious that a small bubble will collapse more rapidly 
than a large bubble. Furthermore, a spherical bubble will collapse more rapidly 
than a bubble of another shape due to a greater surface area-to-volume ratio. Trapped 
vapor will be spherical in shepe for all diameters less than 0. 64 cm (0. 25 ini, the 
channel width dimension. All bubbles of a greater volume will be constrained by the 
channel width and, in the extreme, will acquire the channel shape. Vapor bubble 
condensation time can be determined from Equation 3-61. No comparable equation has been 
identified for determining collapse time for a slab or block of vapor. 


The test plan outlined in the follow ng paragraphs is intended to provide experimental 
data to model the critical elements described above. This te^t plan considered both 
ground based testing and space based testing. For ground based testing, candidates 
for test types were bench tests, drop tower tests and aircraft tests. Spacelab testing 
was considered for space-based test mg. 


Initially, the ground based testing will be described. 
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6. 1 GROUND BASED TESTING 

The ground based testing will consist of filling of a scale model test tank in order to run 
prechill and fill tests. Vapor collapse tests will be run using a transparent dewar. 

Testing of the vent/purge procedure to remove condensible pressurant from the tank 
was considered as a candidate for testing. The decision was not to run tests since the 
recommended purging and venting operations are well understood and do not require 
any additional technology information. The vent pressure, ni mber of purge cycles and 
mechanics of venting and purging can all be determined analytically. 

For tank chilldown, testing is not required because prechill will lower tank temperature 
to the point where excessive pressures will not occur. Tank chilldown will essentially 
be replaced by the prechill period. 

Drop tower and aircraft facilities were evaluated for providing required refilling data. 
These facilities were judged to be inadequate due mainly to the short low gravity test 
period obtainable with either method. Some prechill information could be obtained 
using drop towers however the thermodynamic analysis and modelling techniques 
indicate that normal gravity testing will provide sufficient information to verify the 
analytical models. 

I* 1 INFLOW TESTING. Inflow testing, consisting of prechill and fill testing will 
be conducted in a scale model orbit to orbit shuttle v OOSi tank shown in Figure 
(>-l. This tank is 74 cm ini in diameter with a midsection length of 
cm t9."' ini. The tank mass is 39.5 kg v?7 lb^i .uui the tank volume is n. 29 m : * ^n.3 
fi**i. This tank was selected over a smaller spherical tank because of its lower tank 
mass to volume ratio. The OOS tank has a window in the cylindrical portion of the t:uik 
that can be used for viewing the experiment. 

ti. 1. 1. 1 Prechill Tests. The objective of these tests is to verify the suitability of 
using Equation 3-11 to model tank prechili. Tests will be run with cold GN" and GIF’ 
into the test tank that is initially at room temperature. The prechill test set up is shown 
in Figure t>-_\ The test tank will be insulated with multilayer insulation .11 LD and 
mounted within a 91 cm pit? ini diameter vacuum chamber in order to minimise external 
heat leak during the test period. 

Ptechill o the tank will be used to eliminate the need for venting a receiver tank containing 
liquid and vapor in low gravity. Preehill can be accomplished by introducing liquid or 
vapor into the propelkuit tank. For the application studied the vapor normally generated 
in chilling the transfer lines to the tank will also be employed to chill the tank. 

Consequently, the prechill tests will be run with low temperature vapor. 

It is anticipated that system preehill will be conducted in the following manner: cold fluid 
vill be Introduced through a spray nozzle for a fixed time period, t >nee How is terminated 
heat exchange will occur between the tank wall and fluid until temperature equilibrium 
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occurs. The tank will be vented down to near zero pressure In preparation for a second 
period of vapor inflow. The process will be repeated as many times as required to 
reduce tank temperature to an acceptable level for no-vent chill and fill. 

Ground tests can be conducted to obtain prechill data applicable to a low gravity or zero 
gravity environment only if one-g free convection heat transfer effects become negligible. 
This will be achieved with the oroposed test program by providing a high inflow vapor 
velocity to create forced convection heat transfer as the dominant heat exchange mode 
with the tank wall. (Only a small percentage of the total heat removal occurs during the 
process of liquid evaporation). 

A test matrix for the GH 2 and GN 2 tests is given in Table 6-1. This test series will 
explore the variables of inflow configuration, vapor flowrate, initial tank temperature 
and, of course, fluid properties. The primary inlet configuration, nozzle 41 in Table 
6-1, will be a spray nozzle configuration. The high degree of agitation, mixing end 
liquid-to-vapor heat transfer with a spray nozzle should make it an ideal inlet configura- 
tion for both the prechill and fill periods. Alternately, a bare pipe will be used as 
nozzle 42. Both configurations will also be tested in the fill tests. Predictions indicate 
that the bate pipe may provide the required chilldown efficiency and agitation to permit 
the simpler bare pipe configuration to be used. 

Instrumentation to be used in the proposed testing are listed in Table 6-2. Figure 6-3 
schematically Illustrates the instrumentation locations. Temperature measurements 
will be made with both thermocouples (copper/eonstantan or chromel/constantan) and 
platinum resistance thermometers (Rosemont Model 134L). The platinum thermometers 
will also serve as liquid gas detectors by increasing the excitation current. All 
pressure measurements will be made with full bridge strain gage type pressure 
transducers i$tatham .Model 350). Flow measurements may be made with turbine 


Table 6-1. GH 2 and GNo Prechill Test Matrix 



— 






CONTROLLED 





RUN NUMBER 


1 ! 

0 

i 3 

4 

j 5 . « L 7 

j 3 9 1C 

PARAMETER 

Inlet Configuration 

1 

1 

1 

1 

1 

1 2.2 

! 1 

■ 1 

| 1 11 

Vapor Flowra'e* 


• 

|5<r 
; 

20*' 

" j 3 O'" 30*' ' 20*' 

! 100' 6C“ ! 20“ 

Test Fluid 



Gaseous 

Mvdrogen 

Gaseous Nitrogen 


Test Tank Initial Temperature 294. 4-K *530R). 

Test rank Initial Pressure = 6.9 kN m- (1.0 psia). 

Entering Vapor Temperature - 77. 7K (140R) 

Terminate Test at 413 kN m- ( t o psia) Tank Pressure. 

* ’ 007- flowratt* condition of 0.004 kg sec *0. 009 lb^ sec) for GH-j and 0.024 kg see 
( 0 . 054 lbjn sec) for GNo 
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type volumetric meters (Flow 
Technology Omniflow Model FTM-5- 
LB), orifice and pressure transducer 
combination head meters, and/or hot 
film anemometers (Thermo Systems 
Inc. Model 1053). A magnetic tape 
digital recorder (Dymec Model 150) 
or a digital punch tape recorder 
(Dymec Model 2400) will be the 
primary recorder for initial and final 
steady state data and the slowly 
varying test run data. An analog 
tape recorder (Ampex Model 1200 or 
1300) will be used to record any 
higher frequency transient data 
required during the test runs . 

6. 1. 1. 2 Tank Fill Tests. Propellent 
tank fill should be conducted in such 
a manner that near-thermooynamic 
equilibrium will exist. This repre- 
sents the condition of minimum tank 


Figure 6-3. Instrumentation 


r* 




increased between the phases. 


liquid and vapor will be approached 
as heat and mass exchange is 

• Tbis condition will be promoted by creating a highlv 

■stated, high veloc'ty condition within the tank during fill. The high v.loolty condition 

T, r.T 7° 1 First, a finer spray (l.e. smaller drops) will occur 

iquid velocity is increased. The reduced drop sizes will enhance heat and mass 

exchange with the ullage. Second, an increased liquid velocity will further improve heat 
and macs excaange with the ullage, of course, fill tests conducted on the ground will 
become more representative of filling in space at higher liquid v.l, antics as free convert. 
Ion effects become less significant. 


to order to promote agitation in the tank, two inlets are used to direct flow towards the 
l °? „ ! [ J Calculations indicate that if the upward facing inlet is placed at 3<T r 

° . * be *f ht of the tank below the cank t0 P the spray nozzle should effectively 
agitate the fluid to allow mixing until the tank is essentially full. One of the primarv 

objectives of the testing will be to evaluate the degree of ullage agitation as the liquid 
level is increased above the upper facing jet. The esoerimental data on ullage 
temperature and pressure versus liquid level and spray nozzle outflow rate will be 
correlated with the analyses presented In Section 2< 3. 


Thefill tests are applicable to both cryogenic and noncryogenic refiling. The scale 
model Centaur-OOS tank shown in Figure 6-1 will also be employed for fill testing. 
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for LH 2 except that initial tank temperature is^otTte^tvarilble!’ *“ t0 ^ 

6.1.2 VAPOR COLLAPSE TESTS Van«, . , 

during the propellant fill process. ' once tanTi ? trapped 111 the acquisition device 
introduced to condense this trapped vapor As w a ° a mpleted ' pressura ni *111 be 
trapped within the acquisition device condensaflon will o ° n0 “ 00nde ” 3ibla rapor is 
vapor cannot remain at the elevated pressure while . ° U ^ because the Propellant 

The only uncertainty with this method? th.? SUrrottnded by subcooled liquid. 
Consequently a model was ,„rm„?“ l 2Z\TT ** ^ 
convection effects must be less significant than coii ' t -° r grou ” d tests t0 be successful 
cannot be eliminated from the pmSuill S, " Ce ,ree on 

appended to the prechill and flu test setup s ' vapor c °Uapse tests cannot be 
recommended. Setup - to3 ' aad . the following test configuration is 

A glass dewar as shown in Figure 6-6 win h« j 

interface while minimizing conduction from L ^ VieWlng of the li( ? uid vapor 

to occur at the liquid/vap^r interface. rZ t , ? ^ Heat transfer is flowed 

the vapor trapped within the capillary device.' " ithin the deWai repre ^nts 
because of the relatively flat shaDe nrnHn^ ^ k 1 /vapor contact area is known 

case because there will be negligible cooluL ‘ eS ‘ S :MuId be a worat 

, a screen acqnls.t.on devic/wLe tb^^^^^Pared 

zt™ ^zzvz^: r ~ 

a.r condifloned blcchbcuse ta order to closely control “ 
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Table 6-3 


LHo Fill Test Matrix 


CONTROLLED PARAMETER 




RUN NUMBER 


— 1 

i 

1 

2 j 3 

4 

! s ; 5 

7 1 3 

1 

9 

lu ! 

i 

NOZZLE 

’ 

i i 

! 

. ! . 


! ! 

! 1 i 1 

i 

. 1 . 

j 

o 

t 

2 i 
i 

LH 2 FLOWRATE* 

1 100% 

■ 

90% | 30% 

i 

70% 

i 

, TBD 30% 

t 

80% 400% 

so% 

30% ! 

ENTERING LIQUID VAPR 

: 101. 3 

101. 31101.31101. 3 

l 1 

lj 101.3' 172. 3' 

101.31101.3 ; 

101.3: 

1 7 2 . 3 

PRESSURE, kN/m 2 (PSIA) 
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a4.: 

u 1 14. 71 (25. 0 ) 

j ' 1 
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) i 
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413 . 
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.40) :<60) 

! 
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.60) j 

i 


INITIAL TEST TANK TEMPERATURE IS LH., TEMPERATURE 
TEST TERMINATED AT APPROXIMATELY 90 PERCENT PILL 

* 100% FLOWRATE CONDITION IS . 18 kg/ sec (0.4 lb m ,’sec) 


Table 6-4. Freon MF Fill res. Matrix 
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RUN NUMBER 
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* 10C% FLOWRATE CONDITION IS 1. 13 ke'aec <2 . 5 ’.bm 'sect 


Table 6-3. Proposed Test Matrix for Vapor Collapse Tes.s 


CONI ROLLED PARAMETER 

1 

n — 

1 - 

1 

FI 

— 

4 

1 

i’ 

? 

T 

< 

TEST FLUID 

uu 

1 

i uii 

Lll.» 

UU 

* FREON 

r ri eon 

T 

i FREON 



! 



! MF 

MF 

; m 1 

INTmL . kJI'tD VAPOR PRESSURE. 

lot. a 

i I01..1 

10!, J 

t tT2. J 


101.3 

101.3 I 

KN/m.- p.sUi) 

, 14 . r) 

i (U.T) 
! 

1 14.7) 

j . 95 ; 

j <!■* 7) 1 

1U.7* 

1 1 14. ?> 
1 

TANK i PINCHE.VJE, 

31. 5 

j dH.s) 

' 137. S 

1 

t)9. i) 

L. 


1 127. t) 

KN/m’ tpuWI) 

(5) 

j ! 

(20) 

(10> 

j ■” | 

*10) 

! ,20) 


6-10 




J .y ■«! f ' V" r" J 


J 1 * i — I 


I i 


J 


J 


GHe Source 


Vacuum ■ 
Source 



Description 

Temperature Sensor, Tank Ullage 

Temperature Sensor, Dewar Inner Wall 

Temperature Sensor, Inside Dewar 

Temperature Sensor, Test Fluid 

Capacitor Type Liquid Level Sensor, Inside Dewar 

Capacitor Type Liquid Level Sensor, Outside Dewar 

Pressure Gage (Typical) 

Pressure Transducer, Tank Ullage Pressure 


Figure 6-5. Typical Vapor Collapse Test System 






1* Evaouat8 th « test tank and backfill with GHe to remove moisture and air. 

2. Evaucate and backfill with test fluid vapors to remove GHe. 

3. Fill the tank with test fluid to the initial level. 

4. Allow temperature stabilization. 

5. Pressurize the tank to the desired pressure with GHe. 

Monitor vapor collapse inside glass dewar visually and record liquid level, 
pressure, and temperature data versus time. 

6.2 SPACED BASED TESTING 


One of the mato features of the orbital refilling approach advocated in this study is that 
the principal features of the refilling process can be verified with ground testing. Thus 
the benefit of space based testing would be a combination of; demonstration of the 
complete operation of the system and verifying that the gravity dependence of the 

Pr0< ;! S i e !, lnV n 1Ved d0eS QOt CffeCt ° Vera11 svstem ^^ton. A spacelab experiment 
would be ideally suited for obtaining this information. Several receiver tank shapes 

and sizes should be studied in order to cover the full range of possible orbital 
refilling applications. 


Spacelab refill tests will include the major processes of draining a partiallv full tank 
purging and venting to remove the noncondensible pressurant, prechill or chilldown. ’ 
tilling, and vapor collapse using pressurant. Liquid hydrogen is suggested as the 
test fluid for reasons similar to those given for its selection In Reference 6-1 . The 
test svstem would be designed to fit on a single •*, acelab pallet. Fssential elements 
of the system are shown In Figure 6-6. 


The system will consist of a supply tank with a channel type capillary device for 
acquisition of liquid in low gravity, several receiver tanks with partial (start baskets) 
or total control (channel) type start baskets, and pressurant storage bottles. Plumbing 
Instrumentation and telemetry hardware will also be required. The test will be 
designed as a complete package however, much of the instrumentation such as data 
acquisition and transmitting systems would be usable for other spacelab experiments 

A portion of the experiment package showing only one of the potential receiver tanks ’ 

Is shown in Figure 6-6. 

Instrumentation should Include: 

1. Temperature sensors at selected locations on the test tank walls, in the tank 
Interior, in the inlet line, outlet line, vent line, and screen acquisition device. 
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Figure 6-6. Refill Systems Test on Spacelab 
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3. Devices to monitor liquid position in the tank Interior. 


4. Pressure measurement devices in the test tank, screen acquisition device, inlet 
line, and outlet line. 

5. Acceleration measurement devloes on the test tank in three axis. 

6. Data acquisition system, preferably with capability of transferring data to a 
receiving station on earth for analysis between test runs. 

Inflow testing would consist of chilldown and fill testing for the selected inlet configura- 
tion at several flow rates. Pressure, liquid position and temperatures would be 
measured. If venting is required, vent flow rate and inlet and exit quality from the vent 
system would be monitored. Each of these tests would be followed by a vapor collapse 
test conducted by injecting pressurant into the receiver to condense vapor trapped in the 
start basket. Tests to determine whether helium could be removed from a partially 
full tank would be conducted with tanks less than 10 percent full. After removing the 
residual liquid, purging and venting will be undertaken in order to remove the helium. 

A complete test series would consist, typically, of draining, purging and venting, 
chilling, filling, and condensing trapped vapor by pressurizing. 
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APPENDIX A 








HYPRS COMPUTER PROGRAM If 

The HYPRS computer program was developed with LRAD funds. It is summarized herein L| 
to provide a more complete document of the study results. ij 

The program models the filling of a tank with a liquid propellant which enters as a uniform 
spray. Heat and mass exchange between liquid and vapor and between the “ d 1 

the tank walls, and the resulting temperatures and pressures are calculated. The model | 
can also include a fixed mass of helium pressurant. j| 

The program is diagrammed in Figure A-l. The main program, HYPRS, reads data and | 
inmLzS^riables. DIFE3 calls DERIV, which contains the computation^ routines, and | 
CONTROL, which controls output; it also integrates the time dependent vari ^ le ®’ USiDg | 
values passed from DERIV. In addition, propellant properties are accessed as Block Data.| 

A more detailed flowchart of subroutine DERIV is shown in Figure A-2. Subroutine 
VAPOR initializes mass, pressure, and temperature, and PROP computes vapor and 

liquid properties. The other routines handle some aspect of the system heat an m £ 

transfer. The model can be broken into four major segments, roughly chronological * 
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Figure A-l. Computer Program HYPRS 






Figure A- 2. SUBROUTINE DERIV 
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tLit^a^ incoming^^ € and^hat'all heat for evaporation comes from 

JSTJSs. WhenVage pressure is greater than or equal to incoming U£d njor 
pressure, but no liquid has accumulated, subroutine DROP is called. The routine 
uses incoming liquid mass, velocity and properties, and computes forced conyectioi 
from the ullage to uniform, spherical droplets. When liquid begins to accumulate, 
subroutine FILL is called. This subroutine uses the same equations to m^e vap - 
liquid heat and mass transfer, however liquid properties are assumed to be the result 
of perfect mixing between the incoming spray and the liquid "Siding » * e * 

the volume of liquid which resides as a spray is assumed to be limited to a percentage 
X ™Ze voLe. During these two plumes of the model, two other std.rout.nes, 
WALL and BALANCE, may be called. S the waU temperature ia above Wquid saturation 
temperature, subroutine WALL is called. This routine determines the area of dry 
and wetted wall and computes convective heat transfer from the dry wa area 
IgT, ^Uin“he.t transfer across Urn wetted area. Once the wall 1 ms reached 
liquid saturation temperature, it is assumed that virtually all the heat s ee 
extracted from the wall, and the subroutine is henceforth by-passed. If the s> 
ha^not reached thermodynamic equilibrium, the energy balance is made by camng 
•subroutine BALANCE. This subroutine computes the ullage energy an mas , 
STd and temperature derived from these q^ities with die pressure 

and temperature initially estimated and used throughout the calculations. If the 
"r* not acceptably close, the- pressure and temperature estimates are 
refined, and subroutine PROP is called again, to recompute properties, and start 
the computational cycle over again. The fourth major model segmen . ts 
h, subroutine EQUIL. when the system is judged to have reached thermal euutlibrtum. 
This WO uld correspond to the phase of tank fill during which vapor ts entrained in the 
liquid, and incoming spray keeps the system well mixed. 

This subroutine supersedes all previous ones; it computes va P° i r -‘^ U “ 
taansfer, ullage energy balance, sad checlrn (or 

the four major model segments and subroutines WALL and BALANCE are prese 
in greater detail. 



A. 1 SUBROUTINE LOPRS 


H liquid vapor pressure exceed, ultage pressure during the initial period of liquid 
flow into the* tank, subroutine LOPRS is called. This subroutine assumes that a 
smaU quantity of liquid will instantaneously evaporate as drops pass through 
ullage, and the remaining liquid will be chilled in the process. This »^cooled 
liquid will strike the tank and evaporate as a result of heat exchange with the walls. 
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This process has been simplified by assuming that ali liquid evaporated is a result 
of heat exchange with the tank walls and, the evaporated liquid enters the ullage as 
a saturated vapor. For these assumptions, the rate at which heat is extracted from 
the wall is: 
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wall 
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where: 

Q - rate of heat transfer from the wall 
wall 

m * liquid mass Inflow rate 
L 

h * latent heat of evaporation, at incoming liquid conditions 

Because all liquid is assumed to be evaporated, the rate at which liquid is added to 
the ullage is equal to the liquid inflow rate, that is 


• « 
m = m , 
v *. 
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and the rate at which energy is added to the ullage is then 
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where: 

E .. - rate of energy addition to the ullage 

ull 

h - enthalpv of saturated vapor 
v- sat 


An ullage pressure, P u , is estimated, and the ullage temperature, T u is looked 
up as a function of the estimated pressures and computed energy. P u is then compared 
to P u ’, ullage pressure computed using the equation of state 
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where: 
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R • g*s constant 


m * vapor mass 

v 


V * ullagt? volum* 
a 


„ o,, , w pr .,.ar.. u. no, .ccpubly «1— . U» pro.»ur» aatlmat. » "vtaad and «» 
process re- Iterated until convergence is achieved. 


A. 2 SUBROUTINE DROP 


Wh« ullaga pranaura l* at laaat aa graa, aa liquid vapor pwwra. oalwUtwna 
mada on tha baala of a typinal or avaraga aphartaal droplat. Sibnaititio DRa f 
tnaorporataa Utla modal lor tha oaaa >1 no gaaumulatad liquid. Tha typical . J* 

S. ba uniformly a. tamparatnra axaap, tor a aurf.ca layar o, nag l«Mr 
thicknaaa at taraparatura T„ liquid auturatlon tamparatnra .wraapanding tonUaga 
uraaaura. Furihar, tfta dmplat la aaaumad to ntova at constant valoaity • ' • h 
;; u ” ,m,ou valoclty. Haat tranalar la modnltod aa Inroad corn-action mm, tba u Inga 
to tha droplet, and conduction from tha surface to tha Intarior ot tba droplet. rh. 
convective heat transfer coefficient ts given by 
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where: 


- thermal conductivity 


* average droplet diameter 


incoming liquid velocity 
specific volume of vapor 


viscosity of vapor 


specific heat at constant pressure of vapor 
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Figure A-3. Temperature Distribution within the Droplet. 


The total rate at which heat is convccted to the droplets is then 
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where: 


~ average droplet surface area 

n - number of droplets in the tank during the computation interval. 

The a\erage temperature change within the droplet due to conduction from the surface 
to the interior of the droplet is modeled as a function of the dimensionless parameter 
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where: 

k * thermal conductivity of liquid 
£ 

At =» computation time interval 

v = specific volume of liquid 

c a specific heat of droplet 

r x average droplet radius 

This function, represented here as t\x) , is incorporated into the program as a 

poly nominal cur.e fit to the curves derived in Reference A- 2, and shown in Figure A- 3. 

The average droplet temperature change is given by 

AT * iT s - T d > x f tx) , tA-St 

and the total rate at which heat is conduced within the droplets is: 


«KD ’ m < NC i ' AT 

Evaporation is assumed to occur when exceeds Q ; in this ease the 
. D KD 

evaporation rate, n\^ is: 
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where: 

h. - latent heal of evaporation at ullage pressure 
tg 

h - enthalpy of liquid saturated at ullage pressure 
s 

hp r incoming droplet enthalpy 


If the droplets are not totally evaporated, the change in energy of the remaining liquid 
is given by 





If, on the other hand, is less than or equal to 0^, condensation occurs, and 
the condensation rate, topi is caressed as 
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while the change in droplet energy is 
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A. 3 SUBROUTINE WALL 


Anv liquid not evaporated in the ullage then strikes the tank wall. U the wall 
temperature is greater than liquid saturation temperature, the droplets will flatter 
and boil. If the droplet mass is not sufficient to completely cover the wall, convection 
the walls to the ullage will take place over the unwetted area. This is modelled 
iroutine WALL. The rate of boiling heat transfer, Q w » is given b> the equation 


irom 
in subroutine 




where: T = wall temperature 

q = liquid surface tension 


= specific volume of liquid 


The convective heat transfer coefficient, wall to ullage is 
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where: 


V_ = tank volume 


eff = is defined by equation 3-10, and is assumed to be 0.40 for this study 


resulting in a hea* transfer rate from wall to ullage of 


h x (T - T„) x A 
w v w u ' y 


(A- 161 


where: 


A = area of dry wall 
w 


The total rate at which heat is transferred from the wall, 6_ is the sum of Q and 
§WG * 8111(1 1116 rate of evaporation at the wall, m^.is given as W 


h - h + C , (T_ - T ) 
f s l S D 
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where: 


enthalpy of liquid evaluated at film temperature, 


(T + T )/2 
w s 


By the time the wall temperature reaches liquid saturation temperature, virtually 
all the energy has been extracted from it, so alter that point, the wall calculations 
are bypassed. 


A. 4 SUBROUTINE FILL 


When liquid begins to accumulate in the tank subroutine FILL is oalled 
and the heat and mass transfer between liquid droplets and the ullage is calculated 
in the same way, but the liquid properties are assumed to be the result of perfect 
mixing between incoming liquid, and liquid residing in the tank. Further, the mass of 
droplets is no longer computed as the mass of liquid entering during the computation 
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interval, but rather as a function of the total liquid mass in the and the ullage 
volume. It is assumed that the only significant heat and mass transfer between vapor 
and liquid is that which occurs between the ullage and that liquid which exists as a 
spray. The internal energy of the mixed liquid, u £ , is 

U 1 ' V ( “i*<V u 4 ,/m £ ) (A- 18, 

where: 

u^ * internal energy of liquid before this time step 
m £ 3 total mass of liquid - that which resided in the tank 
at the previous time step, plus the current inflow. 

The mass oi liquid existing as a spray, m , is defined as 
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(A- 19) 


whre: 

V - ullage volume 
n 

p - droplet packing factor- that percentage of tank volume not occupied 

by bulk liquid which can be occupied by liquid spray. For simulations 
run for this study a p of 0.02 was assumed. 

Using liquid properties based on u ,, and the spray mass, m , heat transfer rates, 
liquid evaporation or condensation rates, and the resulting lfquid temperature and 
energy can be computed. If liquid temperature and energy exceed or are equal to 
that of liquid saturated at tank pressure, it is assumed that the tank contents have 
reached thermal equilibrium. If this is not the case, a different path is followed, 
subroutine BA1-ANCE is called and the ullage energy balance is made using the 
mass and energy terms previously calculated. 

A. 5 SUBROUTINE BALANCE 


\Vht?n equilibrium has not been reached the change in ullage energy is expressed as the 
sum of the heat transferred to or from it, the product of the mass condensed or 
evaporated times the appropriate enthalpy, plus the change in ullage volume times 
pressure. That is, 


AE U - IQ 4t ThihAt PA^ 


i 


i 


1 


A- 10 


A- 20 1 





1V-.T. 


a. 


whers h is equal to th« enthalpy of saturated liquid if oonuwuaatlon occurs and the enthalpy 
of saturated vapor if evaporation occurs. The change in ullage volume is given by 
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and the new ullage mass. n» ul . is given by 


m =* m + (m_ + m > At 
Uj u O u 
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where: 


m^ = ullage mass at the previous time step 


Using these quantities a check is made to see if they are in agreement with the estimated 
values of pressure and temperature used in the computations, if the agreement is 
acceptable, end- at- time step values are tabulated and printed, otherwise, the estimates 
of pressure and temperature are refined and the process repeated. 

A, 6 SUBROUTINE EQUIL 

When thermal equilibrium has been reached, calculations of mass and energy balance 
made in subroutine EQl'IL supersede previous calculations. Total mass, 
nip and energy, Up are computed 


ni * m + m 
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where m , the liquid mass includes both liquid inflow, and that mass residing 
in the tana at the beginning of the time step, and uy. the specific energy, is the result 
of perfect mixing. The vapor internal energy, u v , is that corresponding to the 
beginning of the time step. A system temperature is estimated, and pressure, energy, 
and system quality-ami, hence, masses - are computed. If the energy and mass do 
not match U^, and nip the temperature estimate is refined and the process is repeated 
until convergence is achieved. 
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The program user can specify a number of different termination criteria —percentage 
fill* liquid maaa* wall chill-down, or time— and can also specify through input 
the tank siae and material, and the specific propellant of interest. 
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APPENDIX B 


PRE-CHILL OF CONTAINER WITH COLD GASEOl'S HYDROGEN ^OH 3 t 


A schematic of preohtll with OH 3 is shown In Figure B-L Cold hydrogen will flow in at 
constant total enthalpy constant temperature and pressure*. Heat transfer from the 
tank wall to the QNg will oool the tank wall and heat up the GHj, When the wall tempera- 
ture approaches QHg temperature, tank venting will ooour to expel the high vapor. 

Thla till and vent process will be repeated until the tank wall Is chilled to the desired 
temperature. The time rate of change of energy being stored in the GH 3 is equal to the 
time rate of energy being Input during Inflow or taken out during venting, 
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Figure B-l. Prechill Schematic 
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where 

m 3 mass flow rate 

V 3 flow velocity 

g = standard aea level acceleration 

j = mechanical equivalent of heat 

m mass 

h * specific enthalpy 

h - time rate of ohange of h 

u Internal energy 

u time rate of change of u 

q - heating rate 

p - gas pressure 

VOL 3 volume of container 

Tg * container wall temperature 

o * gas density 

A * heat transfer surface 
\v 
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W w - mass of container 


contain* r wall temperature 


time rat* of change of T w 


Npr - Prandtl number 


viscosity 


Cp * oonstant pressure specific h*at of OHg 


C w * sp*olfte heat of container 


acceleration 


X * charact*rUtlc dimension 


thermal conductivity 


Equations B-l through B~S were used to determine the time tv' prechill a container to a 
specific temperature using cold illl*. The simultaneous solution was performed using a 
digital cv'mpute r program. The program was written In a series of subroutines. 


Subroutine ’PHKSVR” Is the principal subroutine to tie called by the main program 
"HYPRfc'S. " Figure R-2 Is a flow chart for “PHESl'U. " 
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